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I. 


INTRODUCTION 


Changing  the  attitude  of  an  object  generally  relies  on 
a  common  technique:  modifying  its  angular  momentum.  The 
mechanisms  available  to  effect  this  change  differ  depending 
on  whether  the  object  resides  in  a  medium  such  as  air  or 
water  or  if  it  exists  in  a  vacuum.  Everyday  experiences 

such  as  the  turning  of  a  car  or  boat  or  the  banking  of  an 

airplane  have  given  many  people  an  intuitive  feel  for  how 
objects  react  on  Earth.  This  intuition  does  not  translate 
to  objects  in  space  because  the  motive  forces  used  on  Earth 
are  not  available:  there  is  nothing  in  the  environment  to 
push  off  of,  hold  on  to  or  generate  lift  from.  One 

solution  that  has  been  developed  is  to  trade  the  angular 

momentum  of  the  spacecraft  with  a  component  mounted  inside 
the  spacecraft.  This  allows  the  total  angular  momentum  of 
the  system  to  nominally  remain  constant,  conservation  of 
momentum,  while  still  allowing  the  spacecraft  to  spin  via 
momentum  transfer.  These  devices  are  generally  referred  to 
as  momentum  exchange  devices  (MEDs) . 

A.  MOMENTUM  EXCHANGE  DEVICES 

Although  there  are  several  different  types  of  MEDs, 
the  majority  of  those  used  today  can  be  divided  into  two 
groups.  The  first  group,  reaction  wheels  (RWs)  consist  of 
a  relatively  large  mass  (a  rotor)  mounted  on  a  shaft  and 
attached  to  a  motor.  The  motor  allows  the  RW  to  change  its 
angular  momentum,  by  changing  the  rotational  speed  of  the 
rotor,  and  thus  the  angular  momentum  of  the  spacecraft  to 
which  the  RW  is  mounted.  The  second  group  of  momentum 


exchange  devices  is  single  gimbal  control  moment  gyroscopes 

1 


(CMGs)  .  The  basic  idea  behind  CMGs  is  to  mount  a  reaction 
wheel  rotating  at  a  nominally  constant  speed  on  a  second 
shaft  that  can  gimbal .  As  the  second  shaft  gimbals,  the 
direction  of  the  momentum  vector  of  the  rotating  wheel 
changes,  thus  changing  the  angular  momentum  vector  of  the 
spacecraft  as  a  whole.  Note  that  variable  speed  CMGs  can 
be  implemented  to  combine  the  characteristics  of  both  RWs 
and  CMGs  into  one  device  [1,  p.  179] . 

Combining  several  MEDs  in  different  orientations 
results  in  an  MED  array  and  allows  control  of  the 
spacecraft  momentum  along  additional  dimensions. 

B.  OBJECTIVES  AND  SCOPE  OF  THESIS 

When  conducting  an  initial  design  of  an  attitude 
control  system  (ACS) ,  attitude  control  devices  are  selected 
based  on  the  mission  of  the  satellite,  gross  design 
characteristics  (spin-stabilized,  tumbling,  three-axis 
stabilized)  ,  mission  life  and  other  big  picture 
specifications.  Using  back-of-the-envelope  calculations, 
initial  sizing  is  performed,  and  then  a  simulation  is 
developed  to  verify  and  iterate  on  the  decisions.  This 
process  continues  until  a  final  design  has  been  established 
that  meets  various  requirements  such  as  mass,  performance 
and  service  life. 

The  goal  of  this  thesis  is  to  develop  equations  and 
algorithms  to  help  automate  the  process  of  analyzing  the 
expected  performance  of  an  arbitrary  MED  array.  This  will 
assist  in  determining  an  initial  component  size  and 
location  more  accurately  than  back-of-the-envelope 
computations,  based  on  required  performance. 


2 


The  developed  algorithms  are  designed  for  ease  of 

implementation  in  a  vector-  or  matrix-based  language  such 
as  MATLAB  and  can  utilize  the  benefits  of  parallel 

processing  when  available. 

The  initial  analyses  will  consider  the  well-studied 
areas  of  angular  momentum  and  torque.  These  provide 

important  characteristics  of  the  array  performance.  By 
using  a  clearly  defined  system  to  describe  the  orientation 
of  MEDs  inside  a  spacecraft,  universal  base  equations  were 
derived  to  eliminate  the  need  for  the  creation  of  one-off 
solutions  for  a  given  array.  These  equations  are  not 
restricted  by  array  size  or  component  orientation. 

One  issue  with  performing  a  direct  comparison  between 
RWs  and  CMGs  comes  from  the  fact  that  any  torque  trajectory 
for  a  system  of  RWs  results  in  a  continuous  trajectory  in 
state-space.  This  relationship  stems  from  the  fact  that 
torque  is  not  state  dependent.  This  is  not  the  case  for 

CMGs,  as  the  torque  available  to  a  CMG  at  any  given  moment 
depends  on  the  current  state  of  the  CMG  array.  To  overcome 
this  problem,  a  state  space  analysis  for  CMG  arrays  is 
developed  as  part  of  this  thesis. 

The  final  objective  of  this  work  is  to  examine  the 
reaction  forces  and  moments  seen  by  the  MEDs  as  a 
spacecraft  maneuvers.  This  particular  area  has  received 
little  to  no  attention  in  the  literature. 

In  general,  momentum  exchange  between  MEDs  and  the 
host  spacecraft  is  discussed  in  precisely  that  term: 
angular  momentum  of  the  system  as  a  whole  is  conserved  by 
exchanging  angular  momentum  between  the  two  bodies.  The 
mechanism  for  this  exchange,  however,  has  a  physical 
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consequence  in  the  form  of  applied  forces  and  moments  on 
the  attachment  points  between  the  spacecraft  and  the  MED. 

When  operating  under  rigid  body  assumptions  and 
modeling  spacecraft  maneuvers,  these  reactions  can  be 
ignored  as  they  are  conservative.  However,  they  represent 
possible  failure  modes  and  thus  must  be  either  considered 
when  designing  the  connections,  such  as  bearings,  otherwise 
the  connections  may  be  significantly  over-engineered  adding 
to  cost  and  mass. 

The  force  analysis  will  also  allow  for  the  comparison 
between  different  control  laws.  This  can  be  useful,  for 
example,  if  it  is  desired  to  extend  the  life  of  components 
by  minimizing  the  peak  forces  and  moments  applied  during 
maneuvers  through  a  redesign  of  the  attitude  control 
system. 

C.  THESIS  LAYOUT 

Chapter  II  presents  background  material  necessary  to 
understand  the  derivations  of  equations  in  matrix  form 
rather  than  vector  form. 

Chapter  III  describes  the  derivation  of  the  equations 
used  to  calculate  the  angular  momentum  space  for  RW  and  CMC 
arrays.  It  also  presents  the  momentum  inversion  algorithm 
that  can  be  used  to  find  the  gimbal  angles  necessary  for  a 
CMC  array  to  orient  the  maximum  net  angular  momentum  vector 
in  a  given  direction. 

Chapter  IV  derives  the  equations  used  to  determine  RW 
array  torque  space.  For  CMGs,  it  shows  the  background  and 
principles  behind  state  space  analysis,  including  the 
determination  of  singularity-free  momentum  space. 
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Chapter  V  steps  through  the  derivation  of  the 
equations  used  in  the  reaction  force  analysis.  This 
includes  determining  both  the  inertial  motion  of  the  MED  as 
well  as  the  applied  forces  and  moments  used  in  the  force 
and  moment  balance  equation. 

Chapter  VI  presents  several  applications  of  the  topics 
discussed  in  previous  chapters.  These  applications  include 
both  notional  array  constructs  and  telemetry  from 
previously  validated  simulations. 

Chapter  VII  gives  the  overall  conclusions  noted  from 
the  research  and  provides  some  areas  that  would  benefit 
from  future  analysis. 
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II .  MATRIX-BASED  SOLUTIONS 


Although  coordinate  transformations  and  other  vector 
mathematics  associated  with  space  systems  engineering  tends 
to  favor  vector  or  dyadic  based  solutions,  the  equivalent 
matrix  form  will  be  used  for  the  derivations  in  this 
thesis.  The  primary  purpose  for  this  decision  is  to  allow 
for  simple  implementation  in  matrix-based  languages  as 
MAT  LAB .  The  bulk  of  the  following  derivation  is  based 
heavily  on  Hughes  [2,  pp .  22-24] . 

Equation  (1)  gives  one  of  the  most  used  relationships 
when  analyzing  spacecraft  dynamics:  the  transport  theorem 
[1,  pg.  12]  .  It  allows  a  vector  expressed  in  one  frame  to 
be  differentiated  with  respect  to  another  frame. 

=  +  (1) 

dr  ’  dr  ’ 

Equation  (1)  is  generally  used  to  find  the  inertial 
derivatives  of  vectors  that  are  more  easily  expressed  in  a 
non-inertial  frame.  The  equivalent  matrix  expression. 
Equation  (2),  involves  transforming  the  vector  into  the 
appropriate  frame,  then  using  the  product  rule  to 
differentiate. 


^v=lR\ 

^v=''k‘‘v+^R‘‘v 


(2) 


is  the  direction  cosine  matrix  (DCM)  that  transforms  the 
vector  “v  from  the  a  frame  to  the  N  frame. 


Although  Equations  (1)  and  (2)  are  mathematically 
equivalent.  Equation  (2)  is  not  written  in  a  generally 
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useful  form  as  the  time  derivative  of  the  DCM  between  two 
frames  is  not  a  physically  measured  quantity.  This  results 
in  additional  computational  overhead  as  that  quantity  must 
be  calculated.  To  overcome  this,  the  DCM  rate  can  be 
rewritten  as  the  product  between  the  DCM  and  a  skew- 
symmetric  matrix  of  the  rotation  rates  between  the  frames. 

First,  consider  the  non-collinear  frames  a  and  N 
fixed  in  inertial  space  and  described  by  the  vectrices^  [2, 
p.  8]  {a}  and  {b\  ,  defined  as  Equation  (3)  . 


'h 

II 

II 

4 

a, 

J 

The  DCM  [R)  between  the  two  frames  has  the  properties 
noted  in  Equation  (4)  . 

';r;r  =  ;r^r  = 

Now,  allow  the  frames  to  rotate  at  different  rates 
with  respect  to  the  inertial  frame.  These  rates  can  be 
represented  by  the  free  vectors  ^^unertmi  ^biinemai  ^ 

relative  rate  between  the  frames  is  simply  the  vector 
difference  between  the  two  as  given  in  Equation  (5) . 


^  The  term  vectrix  (plural:  vectrices)  was  "coined"  by  Hughes  to 
describe  a  convenient  notation  with  which  to  keep  track  of  the 
multitude  of  reference  frames  used  in  spacecraft  attitude  dynamics. 
They  are  particularly  useful  when  converting  vector  equations  to  their 
scalar  equivalents  [1,  pg.  522], 
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alb  a/lnertial  bUnertial 

(O  —  (O  —CO 

^J)la 

CO  =  —CO 


(5: 


Equation  (6)  differentiates  {a}  in  frame  b  with 
respect  to  time  using  the  transport  theorem  and  gives  the 
expected  result. 


( (/{a} 

^  dt 


^  dt  j 


+  co‘^'^x\a 


d  {a} 

y  dt  j 


(6: 


=  0 


Differentiating  |Z)}  in  its  own  frame  with  respect  to 
time  yields  an  important  equality,  shown  in  Equation  (7) : 


Y 

V 


d{b] 

dt 


dt 


dia 


dt 


(7: 


,Mlb 

CO  = 


CO- 


0 

alb 


-CO, 


alb 


CO^ 


alb 
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-CO. 


alb 


CO, 


alb 


-COi 

0 


alb 


Since  |a|  cannot  be  zero,  this  forces  the  parenthetical 
sum  to  add  to  zero  to  satisfy  the  equality  in  Equation  (8) . 


0  =  +  ^Rco‘^ 


:8: 


r  -  =-:r:r^  =:r:r^ 


Returning  to  v,  consider  v  as  a  free  vector  in 
(Equation  ( 9)  )  . 


V  = 


(9: 
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To  express  the  vector  in  either  frame,  use  the  inner 
product  between  the  vectrix  and  the  vector,  given  in 
Equation  (10) . 


v  =  {ay  “v 

v  =  {bY  \ 


(10) 


This  changes  Equation  (2)  to  the  form  found  in 
Equation  (11). 


b 


d^v 

dt 


d{af  ^ 


dt 


v  +  [aY  “v 


d  YY 


dt 


dv  „iu  „ 

—  +  0)  X  V 

dt 


(11) 


Equation  (11)  matches  Equation  (1),  which  is  an 
expected  result.  The  penultimate  step  of  Equation  (12) 
involves  expressing  0}“'^  in  terms  of  frame  a  and  then  using 
matrix  identities  to  rewrite  Equation  (11)  in  the  same  form 
as  Equation  (2 )  . 


CO 


alN 


'  d 

{a} 

alN 

r  a  alN  “1^  f  ■)  1 

V 

dx 

J 

=  CO 

x{a}  =j 

L  ®  J  {«}} 

WT  a  alN  (  1  f  1  f  '\T  a  a/N 

CO  =|<^/  ® 


M"  =M 


(12) 


Einally,  combining  Equation  (10) -(12)  gives  the 
desired  result  in  Equation  (13)  . 

{bY  V  =  + 

%  =  ‘‘v  +  {b]{aY[‘‘co‘‘'Y‘‘v 
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(13) 


By  inspection  with  Equation  (2),  it  is  clear  that 


Equation  (14)  holds. 


N  ry  a  ^^alN 


:r  =  "r‘‘o}- 


(14) 


A  more  detailed  proof  of  this  fact  can  be  found  in 
Likins  [3,  pp .  107-110]. 

Although  performing  second  order  derivatives  on 
equations  with  several  rotations  can  become  unwieldy  based 
on  the  number  of  terms  generated  with  the  product  rule,  the 
relatively  small  number  of  coordinate  frames  to  be  analyzed 
circumvents  this  problem.  An  additional  benefit  is  seen  in 
the  removal  of  the  cross-product  as  a  function.  This 
eliminates  the  need  for  additional  function  calls  within  an 
equation,  improving  both  readability  and  speed.  It  will  be 
necessary,  however,  to  create  the  skew-symmetric  matrices 
manually,  slightly  adding  to  the  length  of  the  code. 
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III.  ANGULAR  MOMENTUM  SPACE 


The  momentum  space  available  to  an  MED  array  provides 
an  indication  of  how  much  momentum  the  array  can  exchange 
with  the  spacecraft  in  any  given  direction,  stemming  from 
the  fact  that  by  the  law  of  conservation  of  angular 
momentum.  Equation  (15),  angular  momentum  must  be  conserved 
in  the  absence  of  external  torques. 

Ktai  =  Kc  +  Ked  =  const.  (15) 


Physically,  this  "maximum  momentum  exchange"  value 
represents  the  maximum  rate  the  array  can  cause  the 
spacecraft  to  spin  (or  the  maximum  spacecraft  rate  that  the 
array  can  stabilize),  shown  by  Equation  (15) .  The  shape  of 
the  MED  array  momentum  space  is  unique  to  each  individual 
array  and  is  based  on  both  the  orientation  and  the  physical 
and  mechanical  properties  of  each  MED. 


^5/C  ^MED  SIqC^SIC  ^MED 


(15) 


A.  REACTION  WHEELS 
1 .  Overview 

Consider  a  perfectly  constructed  RW  rotating  about  its 
z-axis,  as  shown  by  Figure  1.  A  scalar  calculation  of  the 
rotational  momentum  of  the  RW  is  given  by  Equation  (16)  [2, 
pg.  158],  assuming  the  reference  frame  is  body  fixed  to  the 
RW  and  lies  along  the  principal  axes  and  the  RW  is  fixed  in 
inertial  space. 
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Figure  1.  Typical  reaction  wheel,  after  [4] 


(16) 


Generalizing  to  three  dimensions,  while  still  assuming 
the  same  reference  frame  as  above,  yields  Equation  (17), 
while  Equation  (18)  gives  the  most  general  body  fixed  case 
where  the  principal  axes  and  reference  frame  axes  are  not 
collinear  [5,  p.  40] . 
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Applying  the  assumption  that  the  RW  spins  only  about 
the  body  z-axis  =0)  to  Equation  (17)  yields  the  same 

as  Equation  (1),  showing  that  the  scalar  equation  holds 
only  when  the  body  frame  and  principal  axes  are  aligned  (or 
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aligned  close  enough  to  make  the  products  of  inertia  of 
negligible  magnitude) . 

2.  Converting  to  the  Spacecraft  Body  Frame 

Since  Equation  (17)  only  describes  how  the  RW  changes 
angular  momentum  with  respect  to  its  own  reference  frame, 
determining  what  the  angular  momentum  change  means  to  the 
spacecraft  as  a  whole  is  the  next  logical  step.  Completing 
that  step  involves  converting  the  angular  momentum  vector 
from  the  RW  frame  to  the  spacecraft  body  frame.  Figure  2 
provides  a  spatial  representation  of  the  relationship 
between  the  two  frames. 


S/C  Body  to  RW  Frame 


Figure  2 .  Example  Spacecraft  Body  and  RW  coordinate  systems 

Multiplying  successive  orthonormal  rotations  creates 
the  necessary  DCM.  In  this  case,  the  first  rotation  is 
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with  a  magnitude 
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Using  the  assumptions  identified  above  for  the  RW,  the 
final  relationship  between  the  angular  momentum  in  the  two 
frames  is  given  by  Equation  (20)  . 
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The  angular  momentum  in  the  body  frame  for  the  entire 
RW  array  is  simply  the  summation  of  the  angular  momentum  of 
each  individual  wheel  (Equation  (21)). 

%  =Y%  (21) 

army  nv,  '  ' 

/=1 

Equation  (22)  provides  the  vectorized  form  of  Equation 


(21)  . 
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-sin(/i)sin(;0j)  -sin(/„)sin(;5„)1 

'Krray  =  “  COS  (f;  )  Sill  )  •••  -  COS  ( )  Sill  i  (22) 

cos(;5i)  cos(;^„)  J  4 

Given  a  set  k  of  different  combinations  of  RW  wheel 
rates,  an  array  of  angular  momentum  points,  ’’Krmy  '  c:an  be 

created  by  a  matrix  form  of  Equation  (22),  leading  to 
Equation  (23)  . 

-sin(fj)sin(;^j)  -sin(/„)sin(;^J  0  0 

'krray=  “  COS  ( )  Slll  ( )  •••  -  COS  ( /„  )  Sill  )  0  0  i  •••  i  (23) 

cos(;^i)  cos(;^„)  j[o  0 

Use  of  matrix  Equation  (23)  greatly  speeds  up 
execution  time  in  the  MATLAB  by  obviating  the  need  for 
explicit  for  loops.  This  allows  additional  data  points  to 
be  analyzed  in  a  given  amount  of  time,  allowing  for  a  finer 
analysis . 

Eor  the  more  general  case  where  the  wheel  frame  is  not 
aligned  with  the  principal  frame,  a  slight  modification  to 
Equation  (22)  is  required.  Although  the  RW  still  only 
spins  about  the  z^-axis,  the  coupling  from  the  products  of 
inertia,  as  seen  in  Equation  (24),  results  in  angular 
momentum  being  distributed  along  all  axes  in  the  RW  frame. 
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When  converting  to  the  body  frame,  the  entire  DCM  must 
now  be  included  rather  than  just  the  third  column,  giving 
Equations  (25)  and  (26)  . 


"K..  = 


^  R  ^  R  *  •  *  ^  R 


w  1 
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(25: 


rw  rw  rw 
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In  order  to  create  the  equivalent  form  seen  in 
Equation  (22),  an  augmented  DCM  and  column  array  must  be 
created.  It  will  be  of  the  form  seen  in  Equation  (27)  . 
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b  R  bn 

w  'twi  w  ^rw-) 


DR 
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w  1 
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(27: 


Einally,  to  carry  out  multiple  calculations 
simultaneously.  Equation  (27)  can  be  expanded  as  shown  in 
Equation  (28) . 


barray 


b  R  bn  ,  ,  .  bn 
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(28: 


Creating  the  wheel  frame  angular  momentum  matrix  can 
be  done  efficiently  by  generating  it  row  by  row  (as  in 
Equation  (29))  vice  column  by  column. 


18 


(29) 


yz. 


...  Q, 


Note  that  the  augmented  matrices  only  give  the  actual 
momentum  at  a  specific  instant  in  time.  Including  the  time 
dependency  would  additionally  require  fixing  the  RW  frame 
to  the  wheel  and  keeping  track  of  the  angle  of  rotation 
over  time. 


3 .  Angular  Momentvim  Space 

Due  to  the  fact  that  a  maximum  wheel  speed  exists  for 
each  RW  as  a  mechanical  limitation,  there  is  a  limit  to  the 
angular  momentum  magnitude  that  the  array  can  produce  in 
any  given  direction.  Although  methods  exist  for  the 
analytical  determination  of  the  maximum  momentum  envelope 
[4],  Equation  (22)  provides  a  basis  for  the  rapid  creation 
of  a  large  set  of  momentum  points  within  the  momentum  space 
based  on  a  rigorous  sweep  through  all  combinations  of  RW 
speeds.  A  simple  algorithm  can  then  be  used  to  find  the 
maximal  points  for  all  directions  around  a  sphere  to 
determine  the  angular  space  for  further  analysis. 

A  geometrical  approach  to  creating  this  volume  also 
exists  [6],  however  the  point  cloud  approach  presented  here 
allows  additional  flexibility  when  dealing  with  hybrid 
RW/CMG  arrays. 

The  momentum  space  itself  assumes  that  the  RW  array  is 
fixed  to  the  body  which  is,  in  turn,  not  rotationally 
accelerating  in  inertial  space.  If  the  array  is  not  fixed, 
the  DCM  created  in  Equation  (19)  would  change  over  time. 
Although  this  can  be  accounted  for  by  using  the  appropriate 
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time  history  of  the  angles  associated  with  the  DCM,  the 
shape  of  the  momentum  space  itself  would  change  over  time. 
Similarly,  if  the  array  had  a  rotational  acceleration 
component  in  inertial  space  then  this  acceleration  would 
have  to  be  included  in  the  wheel  speed  term,  affecting  the 
momentum  vector  direction  for  each  RW  changing  the  net 
momentum  vector  direction  for  the  RW  array  as  a  whole,  thus 
affecting  the  shape  of  the  momentum  space. 

Although  this  seems  to  indicate  that  the  momentum 
envelope  would,  in  general,  not  apply  since  the  spacecraft 
is  always  accelerating  in  inertial  space,  the  magnitude  of 
spacecraft  rotation  is  much  less  than  the  magnitude  of 
wheel  rotation  so  the  effect  of  the  spacecraft  rotation  on 
the  momentum  space  of  the  RW  array  is,  in  fact,  negligible. 

B.  CONTROL  MOMENT  GYRO  ANGULAR  MOMENTUM 

1 .  Overview 

Although  the  basic  principles  for  determining  the 
angular  momentum  of  a  CMG  matches  those  used  for  a  reaction 
wheel,  the  extra  degree  of  freedom  of  the  CMG  slightly 
complicates  the  equations.  The  rotor  now  rotates  about  the 
local  x-axis,  which  is  the  CMG  body-centered  direction  of 
the  angular  momentum  vector,  shown  in  Figure  3. 
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Gimbal  motor 
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Figure  3.  Typical  CMG,  from  [4] 


The  CMG  gimbal  rotates  about  the  local  z-axis,  and 
this  rotation  must  be  accounted  for  when  converting  the 
angular  momentum  from  the  gimbal  frame  of  the  CMG  to  the 
spacecraft  body  frame.  This  rotation  is  accomplished  by 
performing  a  z-axis  rotation  in  the  opposite  direction  of 
the  gimbal  angle,  S  ,  prior  to  doing  the  first  rotation  from 
Equation  (19) .  The  additional  rotation  is  included  in 
Equation  (30) . 
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cos(-f)  sin(-7)  0  10  0  cos(-(5)  sin(-^)  0 

-sin(-7)  cos(-/)  0  0  cos(-/?)  sin(-y0)  -sin(-^)  cos(-(5)  0 

0  0  10  -sin(-y0)  cos(-/?)  0  0  1 

cos(^)cos(7)-cos(y0)sin((5)sin(7)  -cos(7)sin((5)-cos(y0)cos((^)sin(7)  sin(y0)sin(7) 

=  cos((5)sin(7)  +  cos(y0)cos(7)sin(^)  cos(y0)cos((5)cos(f)-sin((5)sin(f)  -cos(f  )sin(y0) 
sin(y9)sin((5)  cos((^)sin(y9)  cos(y0) 

(30) 

The  final  equation  for  determining  the  orientation  of 
the  angular  momentum  vector  for  a  CMG  in  the  spacecraft 
body  frame  can  thus  be  determined,  as  shown  in  Equation 
(31)  (or  by  the  more  general  form  of  Equation  (32)) . 

X  0  oitqi  r/A’ 

0  A.  0  0,  =‘R  0  (31) 

0  0  /JLlSj  [  0 

V  =  V,  /,  v„  n,  (32) 

.-A  -A  A  JLhJ  Lw„. 

2  .  Angular  Momentvim  Space 

a.  Momentum  Inversion  Algorithm  (MIA) 

Although  a  sweep  analysis  through  all  possible 
combinations  of  gimbal  angles  would  provide  a  similarly 
maximal  result  as  that  performed  for  the  reaction  wheels, 
the  presence  of  singular  surfaces  [7]  presented  an 
opportunity  to  develop  a  more  unique  approach  that  greatly 
improves  the  speed  with  which  the  analysis  is  conducted. 
Essentially,  this  technique  is  an  inversion  of  the  CMG 
array  angular  momentum  equations  typically  seen  in 
literature,  such  as  Equation  (33)  [8,  p.  466]. 
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(33) 


h  =  ZH.(‘^) 

i=l 

-cos(y0)sin((5j)  -cos((52)  cos  (y^)  sin  ((5,)  cos((5’4) 

=  cos((5j)  +  -cos(y0)sin((52)  +  -cos((5’3)  +  cos  (y^)  sin  (<^4) 

sin(y0)sin((5'j)  sin  (y^)  sin  (<^2)  sin  (y^)  sin  (^3)  sin  (y^)  sin  (<54) 

The  basis  for  the  MIA  analysis  lies  in  the  fact 

that  maximizing  the  amount  of  angular  momentum  each 
in(divi(dual  CMG  places  in  the  (desired  direction  will  yield  a 
maximum  result  for  the  array  as  a  whole.  This  allows  the 

overall  problem  to  be  decomposed  into  n  smaller  problems, 

where  the  direction  vector  is  projected  onto  the  rotational 
plane  of  each  CMG  in  the  array  to  determine  the  appropriate 
gimbal  angle  for  that  CMG  to  maximize  the  momentum  it 

places  in  the  desired  direction.  The  momentum  vectors  for 
each  CMG  are  then  added,  resulting  in  a  least  squares 
solution  to  the  inversion  problem.  Figure  4  shows  the 
result  of  this  algorithm  on  a  four  CMG  array. 
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Figure  4.  Schematic  of  the  inversion  algorithm 


For  the  z'*  CMG,  consider  Equation  (30)  ,  rewritten 
in  the  standard  Ax-b  form  as  Equation  (34) .  In  this  case, 

d  is  the  unit  vector  along  the  desired  angular  momentum 
direction,  A  is  the  unknown  gimbal  angle  vector,  and  ^is 
the  remainder  of  the  matrix  from  Equation  (30) . 


cos[/) 

-sin(;r)cos(y0) 

-sin(;K)sin 

cos  (a) 

sin(7') 

cos(;5)cos(7') 

-  cos  (t')  sin 

(^) 

sin  (a) 

0 

sin(/?) 

cos(y0) 

0 

'  '  '  ■  '  ■  ■  "“v"  ■  '  '  '  ■  '  ^  V 

K  A 

One  problem  that  becomes  immediately  apparent  is 

the  fact  that  in  Equation  (34),  there  is  only  one  angle 

being  solved  for.  Using  standard  methods  for  solving  the 

problem  will  not  work  since,  in  the  absence  of  additional 

constraints,  they  will  treat  the  two  terms  as  decoupled  and 
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solve  one  independently  of  the  other.  Analyzing  this 
problem  in  a  vectoral  manner  provides  a  means  to  move 
forward.  Figure  5  represents  the  initial  problem  for  a 
single  CMG. 


MIA  Analysis 


ro"  CMG  Rotation  Plane 

. Initial  CMG  h 

- d 


Figure  5.  Gimbal  rotation  plane  and  desired  direction 

vector  in  the  CMG  frame 


Geometrically,  the  left  hand  side  of  Equation 
(34)  describes  a  plane  in  the  spacecraft  body  frame. 
Figure  5  provides  a  visual  example  of  this,  with  the  circle 


describing  the 

possible 

orientations 

of 

the  CMG 

angular 

momentum  vector 

and  d 

describing 

the 

desired 

angular 

momentum  direction.  If  d  lies  in  the  rotation  plane,  the 
CMG  can  orient  itself  such  that  its  angular  momentum  vector 

lies  parallel  to  the  direction  of  d  .  If  d  does  not  lie  in 
the  rotational  plane,  then  the  ideal  direction  for  the 

angular  momentum  vector  would  be  along  the  projection  of  d 
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onto  the  plane.  To  do  the  projection,  the  normal  {n*)  and 
anti-normal  (n~)  vectors  are  calculated  based  on  the 
momentum  directions  at  two  orthogonal  gimbal  angles,  such 

as  0°  and  90°  (  A(0°)  and  A(90°)  respectively) ,  shown  in 

Equation  (35)  and  diagrammatically  by  Figure  6. 

=K^{Q°)xK^i9Q°) 

(35) 

«  =i:A(90°)xi:A(0°) 


MIA  Analysis 
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n 


Figure  6.  Finding  the  Rotation  Plane  Normal  Vectors 

The  appropriate  normal  or  anti-normal  vector  is 
chosen  such  that  the  relationship  given  in  equation  (36)  is 
satisfied . 

n*d>0  (3  6) 

The  projection  can  now  be  determined  by 
subtracting  the  perpendicular  component  in  Equation  (37) 
and  Figure  7 . 
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(37) 
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MIA  Analysis 
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Figure  7 . 


Projecting  d 


onto  the  Rotation  Plane 


The  least  squares  gimbal  angle  can  thus  be 
calculated  using  Equation  (38),  as  described  in  Figure  8. 


360°  -  cos  '  ^ilA  (0°)*J^^g^.  j  otherwise 


(38) 
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MIA  Analysis 


Figure  8.  Determining  the  Gimbal  Angle  of  the  Projection 

The  resulting  CMG  momentum  vector  described  by 
placing  the  CMG  at  the  calculated  gimbal  angle  provides  the 
least  squares  solution  for  that  particular  CMG.  Orienting 
the  CMG  at  any  other  gimbal  angle  would  increase  the  norm 
of  the  difference  between  the  normalized  CMG  momentum 
vector  and  the  normalized  direction  vector.  Performing  the 
same  algorithm  on  each  of  the  other  CMGs  in  the  array 
provides  the  individual  least  squares  solution  for  those 
CMGs  as  well.  Summing  the  momentum  vectors  provided  by 
those  solutions  gives  a  net  momentum  vector  for  the  array. 

As  seen  in  Figure  4  there  may  be  some  angular 
difference  (or  angle  residual)  between  the  desired 
direction  and  the  least  squares  net  momentum  vector.  It  is 
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squares  gimbal  angle  solution  for  any  given  direction  in 
which  a  solution  exists  (a  solution  will  not  exist,  for 
example,  in  a  direction  that  is  perpendicular  to  the 
rotational  planes  of  all  CMGs) .  In  effect,  the  preceding 
algorithm  simply  provides  an  initial  guess  so  that  a 
numeric  solver  can  find  a  set  of  gimbal  angles  that  give  a 
near  maximum  net  momentum  vector  for  a  CMG  array.  This, 
therefore,  allows  for  the  construction  of  a  quasi-maximal 
surface  by  sweeping  around  a  sphere  rather  than  sweeping 
through  all  possible  gimbal  angle  combinations. 


b.  MIA  Angle  Correction 

Assuming  unit  momentum,  the  momentum  vector  for 
each  CMG  is  calculated  via  Equation  (39) ,  using  the  gimbal 
angle  given  by  Equation  (38) . 


cos(;i^) 

-sin(/)cos(y^) 

-sin(/)sin(y^)1 

COs((5) 

sin(/) 

cos(;^)cos(/) 

1 

O 

o 

sin  ((5) 

0 

sin(;^) 

C0S(;^) 

0 

Since  all  vectors  are  known  at  this  point,  the 
angle  residual  can  be  determined  using  the  dot  product 
definition  in  Equation  (40)  . 


(40) 


If  falls  outside  of  a  given  threshold,  the  dot 

product  itself  can  be  used  to  solve  for  the  gimbal  angles 
that  reduce  the  residual  to  near  zero  based  on  the  identity 
given  in  Equation  (41)  . 

w  II V 

(41) 

u»v  =  1 
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Therefore,  when  the  angle  residual  is  zero. 
Equation  (42)  will  be  true. 

=  0  (42) 

Rapidly  finding  the  solution  to  Equation  (42)  is 
well  within  the  capability  of  modern  numerical  solvers, 
assuming  the  solution  exists. 

If  the  products  of  inertia  are  non-zero  for  the 
CMG  wheel,  the  angular  momentum  vector  will  be  offset  from 
the  spin  axis  of  the  rotor.  This  will  induce  an  error  not 
taken  into  account  by  the  equations  as  developed.  In 
general,  the  products  of  inertia  will  be  orders  of 
magnitude  smaller  than  the  moments  of  inertia,  greatly 
reducing  the  size  of  the  error.  To  correct  for  this  error, 
the  projection  of  the  angular  momentum  onto  the  rotational 
plane  should  be  used  in  place  of  . 

c.  Momentum  Space  Analysis 

A  maximum  momentum  surface  for  the  CMG  array 
based  on  the  MIA  can  now  be  quickly  constructed.  Similar 
to  finding  the  maximum  surface  for  RW  arrays  given  a  point 
cloud,  the  first  step  involves  creating  a  unit  sphere  with 
the  point  density  determining  the  fidelity  of  the  CMG 
surface  construction.  Iteration  through  the  spherical 
envelope  provides  the  direction  vectors  for  the  algorithm. 

The  output  of  the  algorithm  provides  both  a 
maximal  surface  for  all  analyzed  directions,  as  well  as  an 
angle  residual  which  gives  a  measurement  for  how  well  the 
CMG  array  can  point  its  momentum  vector  in  a  given 
direction . 
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C.  SUMMARY 

This  chapter  developed  the  equations  and  algorithms 
necessary  to  create  maximum  momentum  volumes  for  both  RW 

and  CMG  arrays.  These  momentum  volumes  can  be  used  to 

estimate  the  maximum  theoretical  performance  of  an 

arbitrary  array  and  for  more  advanced  techniques  such  as 

the  analysis  of  hybrid  arrays  and  performance  comparison 
between  gimbal  angle-limited  singularity  free  CMG  array 
momentum  volumes. 
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IV.  TORQUE  SPACE 


A.  OVERVIEW 

The  general  equation  of  motion  for  a  spacecraft  is 
given  by  Equations  (43)  and  (44)  [8,  p.  463],  where  J  is 

the  total  moment  of  inertia  of  the  system,  CO  is  the 
spacecraft  body  rate  vector  and  is  the  MED  momentum 

vector  expressed  in  the  spacecraft  body  frame. 

=  (43) 

+  (44) 

Substituting  Equation  (44)  into  Equation  (43)  and 
assuming  a  torque-free  environment  yields  Equation  (45)  [8, 

p.  463] : 

{^J j  +  X  )  ~  ^  •  (45) 

Equation  (45)  can  then  be  rewritten  to  separate  the 
spacecraft  terms  and  the  MED  terms,  shown  in  Equations  (46) 
-(48)  . 

[  J oJ’'^  +  co^'^  X  J co^‘^  +  co’’'^  X  med^^  ^  ^  ^  ^ 

=^h  +’’h  (47) 

'’h  =  '’h  +’’h  (48) 

Recalling  Equations  (16) -(18)  and  applying  them  to 
Equation  (47)  yields  Equation  (49).  Equations  (46) -(49) 
hold  regardless  of  the  reference  frame  they  are  expressed 
in,  so  long  as  all  terms  are  expressed  in  a  consistent 
frame . 
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(49) 


^MED  ^rS^rw  ^cmg^, 

Equation  (49)  allows  Equation  (46)  to  be  further 
generalized  into  Equation  (50) . 


_  ^ 


T  ^,rwin  T  -  rwIN  j  cmg/N  j  -cmglN  I  j  rwin  j 

+l„.fiO  +1„^„C0  +  0)  X\J„..0)  +K^„C0 


')!  =0 

'  -i  MED 


Equation  (50)  can  be  expressed  by  Equation  (51)  by 
assuming  that  =  0  in  the  spacecraft  body  frame.  This  is 
a  reasonable  assumption  as  reaction  wheels  will  be  nearly 
axisymmetric  about  their  spin  axis,  which  is  fixed  in  the 
body  frame.  Note  that  there  is  also  the  implicit 

assumption  that  the  overall  inertia  of  the  spacecraft,  J  , 
is  much  greater  than  Icmg'  will  remain  roughly 

constant.  Since  the  CMC  terms  are  much  more  complex,  they 
will  be  folded  back  into  quantities  that  will  be  developed 
separately  [8,  p.  464]. 


T/^b/N  .  ^h/N  ^  T^^b/Nl  bV ^  ^  1 

yjCO  +CO  xjco  \sic=  Ur«  +  ^cmg\ 

'  X  )]  =  (51) 


B .  REACTION  WHEELS 
1 .  Derivation 

Although  Equation  (51)  provides  a  method  for 
determining  the  inertial  torque  for  the  MED  arrays, 
analyzing  the  torque  capability  in  the  body  frame  (as  was 
done  for  momentum  in  Equation  (23) )  removes  the 
complication  of  having  to  include  inertial  motion. 


34 


Equation  (52)  provides  the  relationship  between  momentum 


and  torque  within  the  local 

wheel  frame. 

'^h='^T 

(52) 

Recalling 

Equation 

(16),  Equation 

(52) 

becomes 

Equation  (53) : 

w  _  w  r  f 

^^rw/N  ,  j  •  nv/A/' n 

j  • 

(53) 

The  RW  is  axisymmetric  about  its  spin  axis,  so  the 
wheel  inertia  will  remain  constant.  Converting  to  the  body 
frame  is  done  in  a  similar  matter  as  done  for  momentum,  as 
given  by  Equation  (54)  . 


w[_nv  JwL_rw  J 

(54) 

Since  the 

wheel  frame 

is  fixed  with 

respect 

to  the 

body  frame,  the 

DCM  remains 

constant,  yielding  the 

final 

torque  equation. 

Equation  (55 

)  ,  for  a  single 

RW. 

(55) 


Assuming  the  body  is  not  accelerating  in  inertial 
space.  Equation  (55)  can  be  reduced  to  Equation  (56) , 
which,  as  expected,  is  the  time  derivative  of  Equation  (20) 


b 


T  = 


0 

0 

Q 

rw 


(56) 


To  determine  the  torque  generated  by  the  entire  array. 
Equations  (21)  and  (23)  can  be  differentiated  and  rewritten 
as  Equations  (57)  and  (58)  . 
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b 


(57) 


-sin(7i)sin(;0i) 

-sin(r„)sin(;5„) 

0  0  ■ 

Q, 

= 

array 

-cos(;'i)sin(;5i)  •• 

•  -cos(r„)sin(;5„) 

0 

0 

cos{P,) 

cos(y5„) 

0 

A  maximal  torque  space  can  then  be  developed  via  the 

same  method  used  to  develop  the  maximal  momentum  space. 

Since  the  geometry  of  the  RW  array  does  not  change 

with  respect  to  the  body,  the  shape  resulting  from  Equation 
(57)  will  be  similar  to  the  shape  developed  by  Equation 

(23)  . 

C.  CMGS 

Although  the  torque-momentum  relationship  given  in 
Equation  (52)  holds  for  CMGs  as  well,  the  resulting 

derivative  is  more  complex  due  to  the  time  variant  gimbal 

angle,  S  .  This  additional  degree  of  freedom  causes  the 
available  torque  to  become  state  dependent.  In  other 
words,  each  individual  CMG  configuration  has  its  own  torque 
space.  Among  other  problems,  this  makes  visualization 

difficult  as  the  torque  currently  available  to  the  CMG 

array  is  not  a  constant  function  of  the  gimbal  angle  and 
will  change  with  time  and  for  each  gimbal  trajectory. 

Traditionally,  this  issue  causes  most  analyses  to  be 
conducted  on  the  momentum  space  of  CMG  arrays,  and 

significant  work  has  been  done  in  that  area  [  9 ;  10 ;  1 1 ;  12  ;  8 , 
pp .  466-470]  .  In  particular,  much  of  the  work  has  focused 

on  singularity  analysis  and  either  developing  control  laws 
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to  avoid  singularities  or  to  create  CMG  arrays  that 
minimize  singularities  [12].  An  alternative  analysis 
involves  looking  at  the  state  space  of  the  array. 

Developing  a  state  space  analysis  for  CMG  arrays  is 
numerically  intensive,  especially  if  a  high  fidelity  is 
desired.  Furthermore,  there  is  an  issue  of  dimensionality. 
By  definition,  a  state  dependent  plot  of  available  torque 
would  describe  the  momentum  output  of  the  array  for  each 
gimbal  angle  and  gimbal  rate.  When  considering  a  four 
wheel  array,  this  implies  there  are  eight  input  dimensions 
that  must  be  mapped  onto  three  output  dimensions. 

Despite  the  problems  associated  with  working  in  state 
space,  there  are  several  benefits  that  can  be  realized. 
For  example,  by  identifying  regions  of  the  space  that  are 
composed  exclusively  of  non-singular  states  it  is  possible 
to  constrain  the  allowable  gimbal  angles  for  the  array  so 
that  the  basic  pseudoinverse  control  law  is  sufficient  to 
maneuver  the  satellite.  Without  constraining  the  gimbal 
angles,  this  approach  would  normally  result  in  the  loss  of 
control  of  the  spacecraft  after  a  short  period  of  time  due 
to  singularities.  Another  option  would  be  to  select  a 
larger  region  by  allowing  only  hyperbolic  singularities  to 
be  included,  as  these  can  be  avoided  by  standard  techniques 
such  as  adding  null  motion  to  the  pseudoinverse  control  law 
[8,  p.  468]  . 

1 .  State  Space  Analysis 

Equation  (59)  represents  one  of  the  key  assumptions 

for  this  analysis.  This  assumption  allows  the  input 

dimensionality  to  be  reduced  by  half,  greatly  shrinking  the 

scope  of  the  problem.  There  is  still  the  issue  of 
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visualizing  data  that  lives  in  however,  where  n  is  the 
number  of  CMGs  in  the  array.  By  constraining  the  output  of 
the  analysis  to  be  a  scalar  rather  than  a  vector,  the  three 
output  dimensions  can  be  further  reduced  to  one.  This  will 
result  in  the  loss  of  some  information;  however  the 
information  may  be  unnecessary,  depending  on  the  particular 
analysis  performed. 

(59) 

2.  Contiguousness 

At  this  point,  it  is  necessary  to  introduce  the 
concept  of  contiguousness.  In  this  thesis,  a  state  space 
analysis  representation  is  considered  contiguous  if  all 
continuous  trajectories  that  can  be  traced  in  the  state 
space  are  physically  realizable.  When  analyzing  a  single 
CMG,  it  would  require  some  effort  to  present  a  state  space 
visualization  that  is  not  contiguous.  When  reducing 
dimensionality  during  a  CMG  array  analysis,  however,  it  is 
very  important  to  keep  this  property  in  mind.  When 
considering  the  problem  at  hand,  developing  gimbal  angle 
limits  to  achieve  a  singularity  free  operational  momentum 
space  for  a  CMG  array,  a  non-contiguous  analysis  may  allow 
for  a  state  within  this  space  that  can  only  physically  be 
reached  by  violating  gimbal  angle  limits. 

For  example,  consider  a  four  CMG  array.  Using  the 
previously  discussed  assumptions,  the  state  space  of  this 
array  can  be  reduced  from  R"  to  R^  .  To  convert  this  down 
to  R\  the  gimbal  angles  for  two  of  the  four  CMGs  can  be 
represented  by  an  index  number  along  the  X  and  Y  axes  of  a 
topographic  plot.  This  state  index  number  refers  to  a 
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specific  state  combination  of  the  two  CMGs  chosen  for  that 
axis.  Assuming  that  adjacent  state  index  numbers  represent 
gimbal  angle  combinations  that  are  physically  next  to  each 
other,  the  resulting  plot  will  give  a  semi-contiguous 
representation  of  the  system,  shown  in  Figure  9  and  Figure 
10.  The  remaining  two  CMGs  in  the  array  would  be 
represented  in  state  index  two  and  have  a  similar  plot  to 
Figure  10. 


State  Space  Analysis  Sample  Plot 


123456789  10 

State  Index  1 

Figure  9.  Four  CMG  semi-adjacent  state  space 
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state  Index  1  States 


Figure  10.  Reduced  dimensional  state  space  analysis  state 

index  array  with  two  CMGs 

This  semi-contiguousness  arises  from  the  loss  of  one 
degree  of  freedom  in  each  axis.  This  implies  that  although 
any  trajectories  drawn  in  this  representation  will  result 
in  valid  physical  trajectories  for  the  CMG  array,  it  does 
not  represent  the  full  set  of  possible  trajectories.  To 
achieve  this,  each  CMG  has  to  be  able  to  change 
independently,  requiring  a  plot  of  greater  than  three 
dimensions,  which  is  very  difficult  to  visualize. 

Aside  from  indicating  the  problem  with  dimensional 
reduction  and  maintaining  full  contiguousness.  Figure  9  and 
Figure  10  provide  a  guide  for  how  to  determine  gimbal  angle 
limits.  Assume  that  the  plot  in  Figure  9  shows  locations 
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that  result  in  singular  (black)  and  non-singular  (white) 
states  in  the  CMG  array  state  space.  Gimbal  angle  limits 
for  singularity  free  space  can  now  be  pulled  from  that  plot 
by  finding  state  index  limits  that  bound  only  white  areas. 
Specifically,  if  the  gimbals  were  limited  to  those 
represented  by  the  gimbal  angles  in  states  one  and  two  for 
state  index  one  and  five  through  eight  for  state  index  two, 
the  resulting  momentum  space  for  the  array  would  be 
singularity  free.  By  Figure  10,  the  first  CMG  in  the 
array,  as  defined  by  the  user,  must  operate  between  zero 
and  one  radians  while  the  second  CMG  must  remain  at  one 
radian.  An  additional  plot  for  state  index  two  would  be 
required  to  determine  the  gimbal  angle  limits  for  the  rest 
of  the  CMGs . 

D.  SUMMARY 

This  chapter  developed  the  equations  necessary  to 
determine  the  torque  envelope  for  a  RW  array  and  showed 
that  the  geometry  of  the  torque  envelope  will  have  the  same 
shape  as  the  momentum  envelope. 

For  CMG  arrays,  the  groundwork  was  laid  to  use  a  state 
space  based  analysis  to  examine  the  possibility  of 
operating  in  singularity  free  momentum  space.  The  concept 
of  contiguousness  was  established  to  ensure  that  the 
results  of  the  analysis  would  allow  for  the  CMGs  to  reach 
every  state  in  the  defined  non-singular  region  without 
having  to  violate  the  established  boundaries.  Subsequent 
chapters  will  look  at  how  the  limited  momentum  envelope 
compares  to  the  theoretical  maximum  developed  in  the 
Chapter  III,  and  test  the  feasibility  of  operating  within 
this  region. 
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V.  ANALYSIS  OF  REACTION  FORCES 


One  important  area  that  has  received  little  attention 
in  the  literature  is  examining  the  reaction  forces  at  the 
attachment  points  of  the  RWs  and  CMGs  throughout  a 
spacecraft  maneuver.  Current  pedagogy  on  this  topic 
focuses  on  the  transfer  of  momentum  between  the  spacecraft 
and  the  MED.  The  vehicle  of  this  transfer,  reaction  forces 
and  moments  on  the  attachment  points  between  the  MED  and 
the  spacecraft,  is  rarely  discussed.  There  are  several 
possible  causes  for  this.  One  reason  may  be  that 
experimental  data  determined  that  standard  attachment 
design  was  robust  enough.  Or  that  other  loading  such  as 
vibrational  loading  greatly  exceeded  the  expected  magnitude 
of  maneuver  loading.  It  is  even  possible  that  the  analysis 
sits  in  a  particularly  bad  convergence  between  two  fields 
such  that  neither  dynamicists  nor  structural  engineers  have 
a  desire  to  study  the  problem.  Regardless  of  the  reason, 
it  is  important  to  understand  what  is  happening  at  the 
attachment  points  in  order  to  appropriately  size  the 
motors,  bearings  and  gimbals  so  that  total  mass  can  be 
minimized  while  maintaining  an  appropriate  safety  factor  to 
protect  against  failure. 

For  the  analysis  presented  in  this  chapter,  all  motion 
between  adjacent  surfaces  is  assumed  to  be  frictionless  and 
all  forces  on  bearing  surfaces  are  assumed  to  be  acting  as 
a  single  resultant  force  at  the  centroid.  The  effect  of 
gravity  is  neglected.  All  final  reactions  are  defined  as 
acting  on  the  component  being  analyzed. 
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A. 


INERTIAL  MOTION 


Solving  for  the  reaction  forces  for  the  MEDs  begins 
with  determining  the  motion  of  each  component  part  within 
inertial  space.  To  calculate  the  motion  specific  to  these 
components,  data  such  as  satellite  orientation,  gimbal 
angles,  etc...  must  be  known.  The  rates  and  accelerations 
associated  with  these  motions  can  be  estimated  using 
numerical  differentiation  techniques,  but  the  fidelity  of 
the  analysis  may  suffer.  The  most  accurate  results  require 
complete  knowledge  of  the  system  during  the  time  period 
analyzed . 

The  inertial  motion  for  both  the  RW  and  the  gimbal 
frame  of  the  CMGs  is  the  same.  To  begin,  the  inertial 
position  vector  for  the  RW  center  of  mass  must  be  expressed 
as  in  Equation  (60)  . 
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r 

rw 


V  T)  .  b 
^S/C  ^rw 


(60) 


As  noted  in  the  introduction,  there  are  two  equivalent 
methods  to  proceed  from  here.  The  first  method,  shown  in 
Equations  (61)  and  (62),  uses  the  transport  theorem  and 
differentiates  twice,  then  uses  the  appropriate  DCMs  to 
convert  all  terms  to  the  correct  frames. 
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To  maintain  all  vectors  in  a  consistent  basis,  all 
body  terms  must  be  expressed  in  the  inertial  frame. 
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Alternatively,  Equation  (60)  can  be  expressed  so  that 
all  quantities  are  expressed  in  the  inertial  frame  and  the 
differentiation  operations  performed  (Equations  (63)  and 


="P  (63) 

'rw  ^  SIC  ^  'rw  V  U  J  ; 

P+  bR  r^+  bR 

an.=  R+ bR  bR  hR  fn.+ bR  (64) 

=^P+^R'’r  +2(^R^r  )+^R‘’r 

Then,  using  the  identities  in  Equation  (65),  Equation 
(64)  can  be  written  in  the  same  form  as  Equation  (62)  as 
shown  in  Equation  (66) . 
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Equation  (66)  gives  the  same  expression  as  Equation 
(62)  .  The  primary  difference  between  the  two  approaches 
lies  in  the  fact  that  the  product-rule  approach,  the  second 
approach,  requires  that  the  vector  be  resolved  in  the 
inertial  frame  prior  to  differentiation  while  the  transport 
theorem  approach  does  not.  This  can  cause  the  product-rule 
approach  to  create  a  large  amount  of  terms  as  additional 
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frames  are  added  to  the  problem,  making  it  fairly  unwieldy 
to  solve  by  hand.  The  primary  benefits  of  the  product-rule 
approach  deal  with  the  simplicity  in  creating  an  algorithm 
to  do  symbolic  differentiation  on  a  computer  as  well  as 
explicitly  maintaining  all  terms  in  a  consistent  reference 
frame.  In  this  case,  there  are  sufficiently  few  reference 
frames  that  both  approaches  are  reasonable  to  work  with. 


If  the  RW  center  of  mass  is  assumed  to  be  fixed  within 
the  structure,  a  valid  assumption  for  a  rigid  spacecraft, 
then  some  of  the  terms  can  be  eliminated  to  give: 
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(67) 


As  mentioned  earlier,  the  CMG  gimbal  frame  has  the 
same  equations  of  motion  as  the  RW  frame,  rewritten  as 
Equation  (  68 )  : 
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(68) 


For  a  CMG,  the  presence  of  an  additional  rotor  frame 
adds  another  term  to  the  initial  position  vector  and 
requires  an  additional  DCM  to  go  from  the  gimbal  frame  to 
body  frame.  Using  a  similar  analysis  as  above  and  applying 
the  rigid  body  assumptions  yields  the  rotor  equations  of 
motion  (Equation  (69)). 
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When  using  the  transport  theorem  to  develop  Equation 
(69),  an  additional  step  is  required  to  determine  the 
angular  rates,  since  it  is  very  unlikely  that  the  inertial 
angular  rate  and  acceleration  of  the  gimbal  frame  is 
measured.  This  quantity  can,  however,  be  determined  from 
the  spacecraft  body  rates  and  the  gimbal  rate  (shown  in 
Equation  (70) ) ,  both  of  which  are  typically  measured 
quantities  . 


0)^'^  =0}”'^ +0)^"^ 


^g/N  ^  ^blN  ^  ^g/b  ^  ^g/N  ^  ^g/b 


•b/N 


\glb 


=  CO  +  (Z»®  +{  CO  +  CO 


■  b/N  g-,glb  ,  ,giN  ,,glb 

=  CO  +C0  +C0  xco 


XCO' 


g/b 


(70) 


B.  REACTION  WHEELS 

1 .  Case  One  -  Cantilever 

Eor  RWs,  two  cases  have  been  considered,  as 
represented  by  Figure  11  and  Figure  12.  For  case  one 
(which  represents  a  cantilever  arrangement) ,  there  is 
exactly  one  reaction  force  and  one  reaction  moment,  leading 
to  Equations  (71)  and  (72)  . 
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Axial 

A 


Figure  11.  Schematic  of  RW  cantilever  arrangement 
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(71) 


(72) 


The  RW  motor,  assumed  to  be  mounted  with  bearing  one, 
is  not  explicitly  included  on  the  left  hand  side  of 
Equation  (72)  .  It  is,  however,  implicitly  included  in  the 
^Inertial  term  US  the  effect  of  that  torque  on  the  system:  . 

If  it  is  desired  to  use  actual  motor  torque  output  instead 
of  measured  or  calculated  wheel  acceleration,  the 
additional  torque  should  be  included.  The  wheel 
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acceleration  due  to  the  motor  torque,  however,  should  then 
not  be  included  in  the  equations  of  motion.  Doing  so  will 
cause  the  terms  to  cancel  out  and  remove  the  motor  torque 
from  the  overall  reactions. 

The  right  hand  side  of  the  force  balance  equation  is 
straightforward  as  the  acceleration  term  was  developed  in 
Equation  (67)  and  all  quantities  in  that  equation  are 
known.  The  right  hand  side  of  the  moment  balance  equation 
requires  some  additional  analysis  to  be  put  into  a  useful 
form. 

Since  the  rotation  of  the  RW  frame  in  inertial  space 
is  not  directly  measured,  it  must  be  calculated  from  known 
quantities  as  shown  by  Equation  (73)  . 

rw/  b/  rw/ 

+(o/b  (73) 


As  noted  in  Section  A,  to  utilize  the  product-rule 
approach  of  differentiation  it  is  important  to  initially 
express  the  vector  in  the  inertial  frame.  After 
differentiation,  the  frame  may  be  changed  as  desired. 


The  angular  momentum  of  the  RW  can  now  be  calculated 
and  differentiated  (Equations  (74)  and  (75))  .  The  moment 
of  inertia  for  the  RW  is  assumed  to  be  constant  about  the 
spin  axis  . 
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The  third  term  of  Equation  (75)  includes  the  matrix 
•  Physically,  this  matrix  is  the  rotation  of  the  body 


with  respect  to  the  RW  frame,  given  in  body  frame 
coordinates.  Since  it  is  a  skew  symmetric  matrix,  there 
are  two  equivalent  approaches:  change  the  basis  of  the 
rotation  vector  prior  to  converting  it  to  a  skew  symmetric 
matrix,  or  use  a  similarity  transformation  to  transform  it. 
Since  the  data  will  be  given  in  the  form  of  a  rotation  rate 
in  the  RW  frame,  the  first  method  should  be  used.  Thus, 
Equation  (76)  now  contains  terms  which  designate  given  or 
easily  calculated  data. 
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Finally,  converting  into  the  RW  frame  yields  Equation 

(77)  . 
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The  motor  torque  is  equal  to  the  third  component  of 

Equation  (77)  [1,  pp .  179-185].  Under  ideal  conditions, 

the  motor  would  resist  any  change  in  wheel  speed  with 

respect  to  the  body  due  to  the  motion  of  the  spacecraft 

around  the  RW  spin  axis  unless  a  torque  input  is  desired. 
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The  inclusion  of  an  additional  motor  torque  term  in  the 
left  hand  side  of  the  torque  balance  equation  would 
(assuming  ideal  data)  negate  the  spin-axis  torque  and 
corrupt  the  output . 

2 .  Case  Two  -  Fixed-Fixed 

The  solution  for  case  two  (given  in  Figure  12) 
requires  some  additional  assumptions  due  to  the 
overdetermined  nature  of  the  fixed-fixed  system.  These 
assumptions,  along  with  solutions  based  on  continuity  and 
compatibility  conditions  allow  for  determining  the  bearing 
reactions . 


Axial 


A 


Figure  12.  Schematic  of  RW  fixed-fixed  arrangement 


For  forces  along  the  axial  direction,  it  is  assumed 
that  the  non-motor  bearing,  bearing  two,  is  allowed  to 
translate  and  thus  cannot  support  a  load.  This  places  all 
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axial  loading  on  bearing  one.  Similarly,  for  torques  in 
the  axial  direction,  bearing  two  is  assumed  frictionless  so 
all  axial  reaction  moments  are  resisted  by  the  motor  at 
bearing  one. 

For  the  transverse  cases,  reaction  forces  and  moments 
for  each  axis  can  be  calculated  via  superposition  using  the 
solved  cases  in  a  reference  book  such  as  Table  8.1  of 
Roark's  Formulas  for  Stress  and  Strain  [13,  pp .  208-223]  . 

To  place  the  problem  in  the  format  of  the  standard 
tables,  d'Alembert's  principle  as  applied  to  rigid  bodies 
is  used  to  rewrite  Equations  (71)  and  (72)  in  the  quasi¬ 
static  form  of  Equations  (78)  and  (79)  [3,  p.  142] . 

Additionally,  all  forces  and  moments  must  be  resolved  into 
the  RW  frame  prior  to  the  analysis. 
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The  final  reactions  are  given  in  Equations  (80) -(85) . 
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C.  CMGS 

In  order  to  simplify  the  analysis  of  a  real  CMG,  such 
as  the  one  shown  in  Figure  13,  a  model  was  constructed  with 
a  rotor  shaft  mounted  directly  onto  a  gimbal  shaft  (Figure 
14) .  The  only  other  assumption  made  in  this  simplification 
involved  the  end  conditions  of  both  shafts.  All  other 
characteristics  such  as  length,  center  of  mass  location, 
mass  and  moment  of  inertia  were  left  as  variable.  Neither 
the  actual  rotor  nor  gimbal  frames  need  be  in  their 
respective  principal  frames.  Instead,  those  body  frames 
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are  defined  with  respect  to  the  geometry  of  the  system.  If 
the  body  frame  and  principal  frame  are  not  aligned  for  the 
rotor,  the  system  is  assumed  to  be  static. 


Spin  Axis 
Reference 


Gimbal 

Axis 


Figure  13. 


Typical  CMG  Construction,  from  [14] 
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Axial 


BRGl 


CMG 

Figure  14.  Simplified  CMG  model 

The  CMG  analysis  is  broken  up  into  two  parts:  analysis 
of  the  rotor  and  analysis  of  the  gimbal .  The  analysis  for 
the  rotor  matches  the  previous  RW  analysis  very  closely. 
The  gimbal,  however,  adds  additional  complexity  due  to  the 
mounting  of  the  rotor  on  the  gimbal  shaft  in  this  model. 
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Both  the  rotor  and  gimbal  are  independently  analyzed 
as  a  cantilever  and  as  a  fixed-fixed  system,  giving  four 
possible  configurations  for  the  CMG  as  a  whole. 

1 .  Rotor 

a.  Cantilever 

Equations  (71)  and  (72)  both  apply  to  the 
cantilever  CMG  rotor  (shown  in  Figure  15) ,  with  the 
inertial  acceleration  already  determined  by  Equation  (69). 
Equation  (86)  gives  those  rewritten  equations. 
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Figure  15. 


◄ - ► 
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Schematic  of  CMG  rotor  cantilever  arrangement 


All  that  remains,  then,  is  to  derive  the 
derivative  of  the  inertial  angular  momentum  vector  to 
create  the  inertial  moment.  The  angular  momentum  vector 
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itself  has  already  been  discussed  and  is  given  in  Equation 
(87)  . 
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Expanding  the  angular  rate  term  and  writing  the 
vector  in  the  inertial  frame  yields  Equation  (88) . 
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Eurther  simplifications  occur  on  the  realization 
that  =  =  result  in  Equation  (89). 
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Differentiating  as  above  and  converting  back  to 
the  rotor  frame  gives  Equation  (90) . 
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b .  Fixed-Fixed 

The  same  force  and  moment  balance  equations  used 
by  the  fixed-fixed  RW  (Equations  (80) -(85))  also  apply  to 
the  rotor  in  Figure  16.  ,  however  the  transverse  directions 

are  now  the  y-z  plane  and  the  rotation  axis  is  now  the  x- 
axis,  so  only  a  slight  modification  of  the  subscripts  is 
required,  which  are  shown  in  Equations  (91) -(95) . 


57 
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Figure  16.  Schematic  of  CMC  rotor  fixed-fixed  arrangement 
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2. 


Gimbal 


a.  Cantilever 

For  the  gimbal.  Figure  17,  an  additional  force 
and  moment  must  be  considered:  the  reaction  of  the  rotor  on 
the  gimbal.  As  these  quantities  are  known  from  the  rotor 
calculations,  they  simply  add  additional  terms  when  solving 
for  the  cantilever  reactions,  given  in  Equation  (95) . 


Figure  17.  Schematic  of  CMC  gimbal  cantilever  arrangement 
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b .  Fixed-Fixed 

To  compute  the  gimbal  reactions  for  the  fixed- 
fixed  case  (Figure  18),  the  rotor  force  and  moment  must  be 
shifted  to  the  gimbal' s  center  of  mass.  Equation  (96),  in 
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order  to  satisfy  the  condition  that  only  a  single  force  and 
moment  are  applied  to  the  beam. 
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Using  the  idea  of  the  combined  force  and  moment, 
the  final  solutions  are  nearly  identical  to  those  presented 
in  the  fixed-fixed  RW  case  and  are  given  in  Equations  (97)- 
(101)  . 
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D.  SUMMARY 

This  chapter  developed  the  equations  of  motion  and 
applied  forces  and  moments  necessary  to  analyze  the 
reaction  forces  of  MEDs  of  various  endpoint  configurations 
over  the  course  of  a  spacecraft  maneuver.  The  derivation 
of  these  equations  focused  on  requiring  only  the  telemetry, 
geometry  and  operating  characteristics  that  would  be 
readily  available  from  either  an  actual  satellite  or 
simulation . 

The  validation  of  the  developed  equations  will  occur 
in  Chapter  VI,  Section  C.  This  validation  will  show  that 
the  developed  equations  behave  properly  under  various 
assumed  condition.  Following  validation,  simulated 
spacecraft  maneuvers  will  be  analyzed,  with  the  focus  on 
showing  the  impact  that  different  control  laws  can  have  on 
the  reaction  forces  during  similar  maneuvers. 
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VI .  CASE  STUDIES 


In  this  chapter,  a  number  of  case  studies  are  used  to 
show  the  results  of  analyzing  standard  sets  of  both  RW  and 
CMG  arrays  using  the  algorithms  developed  in  the  previous 
chapters.  Additionally,  the  performance  of  non-standard 
analyses  such  as  sensitivity  analyses  on  RW  balancing  with 
respect  to  the  array  momentum  envelope  are  included  to  show 
the  flexibility  of  the  algorithms. 

A.  MOMENTUM  ENVELOPE 

1 .  Reaction  Wheels 

For  the  following  RW  analysis,  the  wheel 
characteristics  are  based  off  of  the  Bradford  Engineering 
W18  RW  [15]  .  The  important  characteristics  of  the  RW  are 
as  follows: 

^Rw  =  +  /-  4000 RPM  =  4\9  rad  /  SQC 
^Rw  =  +  /“  55  rad  ! 
m  =  5kg 

■.0022  0  0 
7^^=  0  .0022  0  kg-md 

0  0  .0045 

=  ±  55  (.0045)  =  .2475  N  ■  m 

a.  Nominal  Momentum  Space 

When  an  RW  array  is  composed  of  a  single  RW,  the 
available  momentum  space  forms  a  line  that  lies  along  the 
spin  axis,  as  shown  by  Figure  19. 


63 


One  RW  Momentum  Space 


Figure  19.  Momentum  space  for  a  single  RW 

For  three  RWs  in  an  orthogonal  configuration,  the 
momentum  space  forms  a  box  where  each  of  the  eight  vertices 
is  obtained  by  operating,  in  one  of  eight  combinations, 
each  RW  at  its  fastest  speed  in  either  the  forward  or 
reverse  direction.  Figure  20  shows  an  example  of  this 
shape . 


64 


Three  RW  Momentum  Space 


Figure  20.  Momentum  space  for  an  array  of  three  RWs  in  an 

orthogonal  configuration 

Adding  a  fourth  RW  and  orienting  all  four  RWs  in 
the  standard  54.73°  pyramidal  orientation  [6],  with  y  and 
P  angles  given  by  Table  1,  creates  a  more  complex  shape:  a 
dodecahedron . 


RW 

7  (deg) 

P  (deg) 

1 

0 

54.73 

2 

90 

54.73 

3 

180 

54.73 

4 

270 

54.73 

Table  1.  Parameters  for  a  four  RW  54.73°  pyramidal 

configuration 
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As  before,  the  vertices  represent  different 
combinations  of  wheels  operating  at  their  maximum  speed  in 
either  the  forward  or  reverse  direction.  The  dodecahedron 
is  shown  in  Figure  21. 


Four  RW  Momentum  Space 


Figure  21.  Momentum  space  for  an  array  of  four  RWs  in  a 

54.73°  pyramidal  configuration 

Although  there  is  more  than  one  configuration  of 
RW  arrays,  the  characteristic  of  a  given  array  will 
ultimately  dictate  its  suitability  for  a  given  application. 
For  the  three  RW  box  array,  each  dimension  is  driven  by  a 
single  RW,  so  the  size  of  each  side  is  equal  to  twice  the 
momentum  capability  of  each  wheel.  For  the  four  RW 
pyramidal  array,  the  skewed  RWs  allow  additional  RWs  to 
affect  the  array  momentum,  even  if  the  desired  direction  is 
along  the  rotation  axis  of  one  of  the  RWs.  Ultimately, 


this  means  the  largest  sphere  that  can  fit  inside  the 

momentum  volume  for  the  pyramidal  array  is  larger  than  that 
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for  the  box  array.  The  theoretical  maximum  performing 
array  according  to  the  maximum  sphere  metric  would  be  one 
that  has  a  perfectly  spherical  momentum  space.  This  would, 
however,  be  impractical  to  build  as  it  would  require  a 
large  number  of  devices. 


Jb.  Skewed  Momentum  Space 

One  benefit  to  calculating  the  momentum  space 
using  Equation  (28)  instead  of  Equation  (22)  is  that  wheel 
imbalances  can  now  be  taken  into  account. 


Assuming  that  the  same  RW  as  analyzed  above  now 
contains  a  large  imbalance  such  that  the  body  z-axis  of  the 
wheel  is  offset  from  the  principal  z-axis  by  10°,  the 
single  RW  momentum  space  line  becomes  skewed  as  shown  by 
Eigure  22  when  frozen  at  a  particular  instant  in  time.  As 
the  wheel  rotates,  the  momentum  vector  will  create  a  cone, 
an  effect  commonly  called  nutation  [5,  pg.53] .  This  makes 
the  momentum  vector  time  dependent  and  complicates  the 
analysis . 

One  possible  skewed  inertia  matrix  is 


2.235x10^'  -3.468x10^' 

-3.468x10  '  2.235x10  ' 

2.781x10^  -2.781x10" 


2.781x10  " 


-2.781x10^ 
4.431x10  ' 


kg-m^  . 
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One  RW  Momentum  Space  (Nominal)  One  RW  Momentum  Space  (Skewed) 


Figure  22.  Balanced  (left)  and  unbalanced  (right)  RW 

comparison 


When  applied  to  a  three  or  four  wheel  array, 
these  unbalanced  wheels  skew  the  entire  momentum  envelope. 
This  is  clearly  shown  in  Figure  23.  These  shapes  only 
apply  when  the  noted  requirement  for  the  applicability  of 
Equation  (28)  is  met:  that  the  shape  is  valid  only  for  a 
particular  instant  in  time.  If  examined  over  time,  the 
shape  would  "breath"  as  each  vertex  moved  with  the  momentum 
cone  for  each  RW  causing  the  faces  to  change  size,  shape 
and  orientation. 

Although  the  10°  offset  makes  the  differences  in 
the  envelopes  clear,  the  offset  represents  an  imbalance 
that  would  result  in  a  shortened  component  lifetime.  More 
realistically,  the  offset  would  likely  be  much  less  than 
1°  . 
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Three  RW  Momentum  Space  (Nominal) 


Three  RW  Momentum  Space  (Skewed) 


Four  RW  Momentum  Space  (Nominal)  Four  RW  Momentum  Space  (Skewed) 


Figure  23.  Balanced  (left)  and  unbalanced  (right)  multi-RW 

array  comparison 

One  question  that  arises  from  this  study  is  what 
effect  this  skewed  shape  has  on  the  total  enclosed  volume. 
A  smaller  enclosed  volume  represents  a  performance  loss,  so 
determining  the  degree  of  this  loss  in  the  range  of 
expected  imbalance  represents  a  useful  exercise.  By 
iterating  the  offset  angle  and  calculating  the  resulting 
momentum  envelope  volume,  a  relationship  can  be 

established,  as  shown  in  Figure  24  and  Figure  25. 
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Three  RW  Array  Imbalance  Analysis 


Figure  24.  Three  RW  array  imbalance  sensitivity  analysis 


Four  RW  Array  Imbalance  Analysis 


Figure  25.  Four  RW  array  imbalance  sensitivity  analysis 

The  results  show  less  than  a  .04%  performance 
drop  at  1°  skew.  The  small  reduction  in  performance  for 
what  is  still  a  large  imbalance  indicates  that  assuming  the 
RW  rotates  about  a  principal  axis  will  have  a  negligible 
impact  on  preliminary  ACS  design. 
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2. 


CMG  Momentvim  Space 


Consider  the  RW  analyzed  in  the  previous  section  now 
mounted  on  a  gimbal .  Further,  assume  that  the  wheel  only 
rotates  at  its  maximum  speed  in  the  forward  direction  and 
that  the  axis  of  rotation  is  the  rotor  x-axis,  to  be 
consistent  with  the  system  developed  earlier.  The 

available  momentum  space  for  a  single  CMG  is  a  circle,  as 
shown  in  Figure  27. 

Figure  28  shows  the  angle  residuals  between  the 
desired  momentum  vector  and  the  momentum  vector  that  the 
CMG  can  produce  via  the  MIA  process  described  in  Chapter 
III,  Section  III.B.2.  The  desired  direction  vector  is 
determined  from  the  azimuth  and  elevation  angles  around  the 
unit  sphere.  Figure  26  provides  the  convention  to 

determine  the  direction  of  d  using  the  (j)  (elevation)  and  6 
(azimuth)  angles.  It  is  clear  that  this  single  CMG  array 
can  only  match  the  desired  vector  when  it  lies  in  the  plane 
of  rotation.  In  other  words,  as  the  elevation  component  of 

d  deviates  from  0°,  the  residual  from  the  MIA  produced 
momentum  vector  will  be  equal  to  that  deviation. 
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4)  and  0  Convention 


Figure  26.  Convention  for  using  (j)  and  6  to  describe  d 


One  CMC  Momentum  Space 


Figure  27.  Momentum  space  for  an  array  of  one  CMC 
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Angle  Residuals 


e(deg)  °  °  ^(deg) 


Figure  28.  Angle  residuals  between  the  desired  momentum 
vector  and  the  actual  momentum  vector 

Placing  three  CMGs  in  an  orthogonal  configuration 
gives  a  much  more  interesting  shape,  seen  in  Figure  29. 
Note  that  this  figure  represents  the  maximum  momentum 
theoretically  achievable  in  a  given  direction,  however  due 
to  the  presence  of  internal  singularities  it  may  not  be 
physically  possible  to  reach  the  surface  in  a  given 
maneuver . 
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Three  CMG  Momentum  Space 


Three  CMG  Momentum  Space 
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Three  CMG  Momentum  Space 
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Figure  29. 


Three  CMG  array  maximum 


momentum  space . 


Although  a  rigorous  singularity  analysis  was  not 
performed,  it  is  likely  that  the  conical  shapes  match  the 
shape  of  a  singularity  surface  [4]  due  to  the  relatively 
large  angle  residual  shown  in  the  vicinity  of  each  by 
Figure  30. 


Angle  Residuals  (deg) 


■!>  (deg) 


Figure  30.  Three  CMG  array  angle  residual  plot 


For  the  54.73°  CMG  pyramid,  the  standard  shape  appears 
in  Figure  31.  The  large  vertical  stripes  present  at  ^=90° 
and  (^  =  270°  in  Figure  32  are  likely  due  to  a  singularity 


present  when 


±1 

0 

0 


In  these  directions,  changing  0  has 


no  effect  on  the  direction  of  d  .  This  means  that  as  the 
MIA  algorithm  sweeps  through  the  azimuth  values,  the  same 
residual  will  be  seen.  Further  analysis  is  required  to 
determine  why  the  vertical  stripe  is  not  preset  at  ^  =  90° 
for  Figure  30. 
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Figure  31. 


Four  CMG  array  momentum  space 
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Figure  32.  Four  CMG  array  angle  residuals 


3 .  Hybrid  RW/ CMG  Array 

By  maintaining  all  momentum  data  in  point  clouds,  it 
becomes  possible  to  analyze  the  maximum  momentum  surface  of 
a  composite  array  with  both  RWs  and  CMGs .  To  the  best  of 
the  author's  knowledge,  such  an  analysis  has  not  been 
reported  in  the  literature.  Figure  33  shows  such  an  array, 
with  three  CMGs  and  three  RWs  oriented  in  two  sets  of 
orthogonal  box  arrays,  similar  to  the  three  device  cases 
analyzed  previously. .  There  is  some  distortion  of  the  CMG 
shape  (beyond  the  increased  size  from  adding  the  volumes) . 
Internally,  the  CMG  singularity  surfaces  get  pushed  to  the 
border  of  the  original  RW  surface,  as  that  provides  a 
singularity  free  volume. 
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This  singularity  free  space  implies  that  it  may  be 


possible  to  include  a  separate  RW  array  on  a  spacecraft  for 
the  purpose  of  elliptic  singularity  avoidance.  The  RWs 
would  allow  a  torque-generating  CMG  motion  to  bypass  the 
singularity  threshold,  potentially  allowing  for  an  expanded 
utilization  of  the  CMG  performance  envelope,  or  the  ability 
to  use  a  simplified  steering  law. 


3  CMG/3  RW  Combined  Momentum  Space 


3  CMG/3  RW  Combined  Momentum  Space  3  CMG/3  RW  Combined  Momentum  Space 


0  5 

(N-m-sec) 


(N-m-sec) 


Figure  33.  Composite  array  momentum  volume  for  a  three 

RW/three  CMG  hybrid  MED  system 
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B .  TORQUE  ENVELOPES 


1 .  Reaction  Wheels 

With  the  momentum  vector  for  the  RW  being  essentially 
fixed  along  the  rotation  axis  (assuming  a  perfectly 
balanced  wheel) ,  the  torque  space  assumes  the  same  shape  as 
the  momentum  space,  as  seen  in  the  torque  space  for  the 
four  RW  pyramidal  array  in  Figure  34.  The  other 

configurations  are  similar  to  the  example  shown  in  the 
previous  section. 


Four  RW  Torque  Space 


Figure  34.  Four  RW  pyramid  array  maximum  torque  space 

2 .  CMG  Torque  Space  Analysis 

a.  Reduced  Dimensional  State  Space  Analysis 

When  dealing  with  state  space  representations  of 
a  CMG  array,  as  discussed  in  Chapter  III  Section  C,  it  is 
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very  important  to  maintain  the  contiguous  or  semi- 
contiguous  condition  when  reducing  dimensionality. 

There  are  several  possible  metrics  that  can  be 
calculated  for  the  state  space  analysis.  One  typical 
analysis  is  based  on  testing  the  singularity  condition  (SC) 
of  a  particular  state  as  defined  by  Equation  (102)  [8,  pg. 

467]  . 

5'C  =  det(^^^)  (102) 

Where  A  can  be  defined  as  in  Equation  (103) . 


Eor  the  case  n>3 ,  a  non-zero  SC  implies  that  the 

pseudoinverse  control  allocation,  A^=A^(^AA^'j  ,  exists 

because  AA^  is  invertible.  Thus,  a  vector  of  gimbal  rates 
can  be  found  that  will  satisfy  Equation  (103) .  Therefore, 
developing  a  state  space  map  based  on  the  SC  of  each  state 
allows  a  region  of  gimbal  angles  for  each  CMC  to  be  created 
that  is  singularity  free. 

Eigure  35  shows  a  map  of  the  SC  based  on  state 
index  value,  with  Figure  36  showing  the  gimbal  angles  based 
on  state  index  number  with  zero  representing  the  bottom 
left  corner  of  the  singularity  free  box  and  increasing  to 
the  right  and  up  for  state  index  one  and  two,  respectively. 
The  state  indexes  themselves  were  created  by  taking  pairs 
of  CMGs  (one  and  two  for  the  first  state  index,  three  and 
four  for  the  second)  and  allowing  the  gimbal  angle  of  one 
CMC  in  each  pair  to  cycle  through  an  entire  circle  before 
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incrementing  the  other  CMG  in  the  pair.  This  continued 
until  the  second  CMG  gimbal  angle  had  gone  through  360°  . 
Figure  37  describes  the  process  in  pseudo  code. 


Raw  SC  for  4  CMG  Pyramid 
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Figure  35. 


Four  CMG 
space 


pyramid  reduced  dimensional  state 
analysis  showing  raw  SC  data 
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Allowable  Gimbal  Angles 


Figure  36.  Four  CMG  pyramid  reduced  dimensional  state 

analysis  gimbal  angle  constraints 


forgimbal2  =  1  to  360 
forgimball  =  1  to  360 
stateindexl(kk)  =  [gimball,gimbal2] 
gimball  =  gimball  +  increment 
kk  =  kk  +  1 
end 

gimbal2  =  gimbal2  +  increment 
end 

Figure  37.  Pseudo  code  representation  of  the  creation 

of  the  state  indices 


Figure  38  shows  the  state  conditions  that  meet 
the  inequality  SC>0.  The  black  box  present  in  Figure  36 
and  Figure  38  indicate  the  identified  singularity  free 
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region,  based  on  a  simple  search  algorithm  that  finds  the 
largest  square  area  that  meets  the  required  condition. 
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Figure  38.  Four  CMG  pyramid  reduced  dimensional  state 

analysis  showing  analyzed  SC  data 

Although  the  gimbal  limits  can  be  pulled  off  of 
Figure  37,  understanding  that  the  steps  seen  in  CMG  four 
represent  the  resolution  of  gimbal  angles  analyzed 
underscores  the  fact  that  this  representation  is  not  very 
accurate.  The  purpose  of  using  such  large  steps  initially 
is  to  reduce  the  analyzed  state  space  since  each  gimbal  has 
to  be  varied  through  the  entire  360°.  By  using  the 
resulting  gimbal  limits  in  a  subsequent  iteration  and 
maintaining  the  number  of  steps  for  each  angle  range,  the 
results  can  be  refined  at  the  same  speed  as  the  initial 
examination.  This  can  be  continued  until  the  desired 
accuracy  is  developed.  Figure  39  shows  the  results  for  the 
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[SC  >  0]  for  4  CMG  Pyramid 
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second  and  third  iteration  while  Table  2  gives  the  gimbal 
angle  limits  for  each  iteration. 


Iteration 

CMG  1 
(deg) 

CMG  2 
(deg) 

CMG  3 
(deg) 

CMG  4 
(deg) 

1 

Low 

0.00 

77 .14 

0.00 

77 . 14 

High 

360.00 

180.00 

360.00 

180.00 

2 

Low 

0.00 

57.30 

0.00 

155.52 

High 

360.00 

118 . 68 

360.00 

229.18 

3 

Low 

0.00 

65 .48 

0.00 

139.15 

High 

360.00 

122 .78 

360.00 

196.44 

Table  2.  Gimbal  limits  per  iteration  for  four  CMG  pyramid 
reduced  dimensional  state  space  analysis  based  on 

SC 


The  point-cloud  momentum  space  resulting  from  the 
calculated  gimbal  angle  limits  is  given  in  Figure  40.  For 
these  gimbal  angles,  there  is  a  much  larger  proportion  of 
volume  in  the  positive  x-y-z  octant.  So  although  the 
momentum  space  volume  may  be  sufficient,  the  location  of 
this  volume  with  respect  to  the  origin  may  not.  For 
comparison  purposes.  Figure  41  shows  the  limited  momentum 
space  plotted  within  the  unlimited  momentum  space. 
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Iteration  2 


[SC  >  0]  for  4  CMC  Pyramid 
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Iteration  3 


[SC  >  0]  for  4  CMG  Pyramid 
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Figure  39. 


Results  of  iteration  two  and  three  for  the  four 
CMG  pyramid  reduced  dimensional  state  space 
analysis  based  on  SC 
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Gimbal  Limited  Momentum  Space 
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Figure  40. 


Momentum  space  for  gimbal  angle  limited  four  CMG 
pyramid  reduced  dimensional  state  space  analysis 

based  on  SC 


A  useful,  though  not  exhaustive,  test  to 
determine  if  the  volume  is  singularity  free  is  to  pick  two 
points  within  the  volume  and  use  a  simple  pseudoinverse 
controller  to  drive  the  four  CMG  array  from  one  momentum 
state  to  another.  Under  normal  circumstances,  the 
pseudoinverse  control  law  tends  to  drive  the  array  towards 
a  singularity,  but  by  staying  within  the  declared 
boundaries  singularities  can  be  avoided. 
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Figure  42  and  Figure  43  show  the  result  of  this 
test,  using  a  simple  PD  controller  to  drive  the  momentum 
error  to  zero. 


Gimbal  Limited  Momentum  Space  Comparison 


h  (N-m-sec) 


h  (N-rti-sec) 


Gimbal  Limited  Momentum  Space  Comparison  Gimbal  Limited  Momentum  Space  Comparison 


Figure  41.  Reduced  dimensional  limited-  and  unlimited- 
gimbal  angle  four  CMG  pyramid  momentum 
space  comparison 


Normal  attitude  control  operations  do  not 
calculate  actual  angular  momentum  values;  they  are  usually 
based  on  current  and  final  quaternions.  If  doing  a 
standard  eigenaxis  maneuver  using  quaternions,  the 
intermediate  step  of  calculating  angular  momentum  and 
ensuring  the  trajectory  stays  within  the  bounds  (difficult 
due  to  the  irregular  shape  of  the  volume)  can  be  neglected 
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by  simply  ensuring  all  CMGs  stay  within  their  gimbal  angle 
bounds . 


Momentum  Error  vs  Time 


Torque  vs  Time 
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Figure  42.  Pseudo-inverse  control  law  operation  in 
singularity  free  CMG  momentum  space 


When  using  the  pseudo-inverse  control  law  as  a 
basis,  there  are  several  difficulties  to  overcome  when 
designing  this  new  implementation.  First  of  all, 
controlling  gimbal  rate  behavior  near  gimbal  angle 
boundaries  is  very  important.  It  would  be  very  undesirable 
to  generate  a  high  gimbal  rate  in  an  adverse  direction, 
forcing  the  CMG  across  the  operational  boundary.  To 
minimize  rate  changes,  one  can  simply  convert  to  a  delta 
form  of  the  control  law  as  given  in  Equation  (104)  . 
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(104) 


5  =  AS  +  Sq 
h  =  a{^AS  +  Sq^ 
AS  =  A^(h-AS,) 
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Figure  43. 


Pseudo-inverse  control 
singularity  free 


law  momentum  trajectory  in 
CMG  momentum  space 


The  properties  of  the  least  squares  solution  will 
then  minimize  ||a^||  ,  thus  minimizing  gimbal  rate  changes 

[16]  . 

It  is  also  necessary  to  ensure  that  the  gimbal 
rates  tend  to  avoid  the  saturation  conditions.  To  make 
this  happen,  a  weighting  factor,  ,  can  be  added  to 
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Equation  (94) .  This  factor  will  act  to  drive  the  solved  AS 
such  that  A5  +  5^  will  direct  the  gimbal  away  from  the  limit. 
The  full  steering  control  law  is  given  in  Equation  (105) . 


h  =  a{^AS  +  + 

AS  =  A^(h-AS  +  Awj.) 


(105) 


Although  no  analysis  has  been  completed  to 
determine  create  a  working  example  of  Equation  (105),  it  is 
likely  that  the  design  of  a  wy  function  to  achieve  the 
desired  behavior  will  be  the  most  challenging  part. 


b.  Full  Dimensional  State  Space  Analysis 

One  of  the  primary  benefits  of  reduced 
dimensional  analysis  lies  with  the  fact  that  intermediate 
results  can  be  represented  in  a  form  that  is  fully 
representable  in  .  A  close  inspection  of  the  state  index 
plot  in  Eigure  39  shows  that  not  all  permutations  of  CMG 
gimbal  angle  states  are  analyzed.  Though  the  states  that 
are  analyzed  are  contiguous,  there  are  additional 
contiguous  states  that  are  ignored.  This  is  the  heart  of 
why  the  previous  analysis  was  characterized  as  semi- 
contiguous.  The  unanalyzed  states  impose  unnecessarily 
restrictive  constraints  on  the  CMG  array  based  on  the 
gimbal  angle  limits. 

A  fully  contiguous  analysis  representation 
requires  a  five  dimensional  plot  for  the  system  being 
analyzed.  This  plot  cannot  be  easily  displayed  for 
visualization  with  current  technology.  fortunately,  the 
same  analytical  principles  still  apply. 
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By  storing  the  SC  values  in  a  four  dimensional 

matrix,  a  hypercube  can  be  formed  within  that  matrix  to 

determine  the  gimbal  angle  limits  for  singularity  free 
space.  With  these  limits,  similar  analyses  such  as  the 

momentum  space  point  cloud  can  be  performed. 

Table  3  provides  the  results  for  each  iteration 
of  the  full  dimensional  state  space  analysis  while  Figure 
44  and  Figure  45  give  the  momentum  space  and  provide  a 
comparison  to  the  prior  methods,  respectively. 


Iteration 

CMG  1 
(deg) 

CMG  2 
(deg) 

CMG  3 
(deg) 

CMG  4 
(deg) 

1 

Low 

25 .71 

102.86 

25 .71 

102.86 

High 

180.00 

257 .14 

180.00 

257 .14 

2 

Low 

0.00 

193.40 

0.00 

193 .40 

High 

113 . 65 

322 .48 

113.65 

322 .48 

3 

Low 

0.00 

167.90 

21.51 

199.92 

High 

86.05 

295.97 

107 . 57 

327 . 98 

Table  3.  Iterative  results  for  the  full  dimensional  state 

space  analysis  of  the  four  CMC  pyramid 

The  complexities  of  the  gimbal  limited  shapes 

make  comparison  between  the  two  methods  difficult.  One 
fairly  simple  metric  to  compare  the  full  and  reduced 

dimensional  shapes  is  to  compare  the  radii  of  the  largest 
sphere  centered  at  the  origin  that  can  be  inscribed  within 
the  volume  of  the  shapes.  For  the  reduced  dimensional 
analysis  case,  this  radius  is  approximately  0.1241/?^^^  N-m-s 
while  for  the  full  dimensional  analysis  case  it  is 

approximately  0.0689^^^  N-m-s.  The  full  dimensional 

momentum  space,  however,  is  much  more  centered  than  the 
reduced  dimensional  analysis  momentum  space.  The  spheres 
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are  shown  superimposed  on  the  respective  shapes  in  Figure 
46  and  Figure  47. 


Gimbal  Limited  Momentum  Space 


(N-m-sec) 


1  0  1 
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Gimbal  Limited  Momentum  Space 


Gimbal  Limited  Momentum  Space 


Figure 


44.  Four  CMG  pyramid  full  dimensional  state  space 
analysis-based  singularity  free  momentum  space 


Although  Crews  ([4])  indicates  that  the  maximum 
singularity-free  sphere  has  a  radius  of  1.0069/z^^g  N-m-s,  a 
number  which  greatly  exceeds  the  normalized  radius  of 
0.0689/z^^g  N-m-s  for  the  full  dimensional  state  space 

analysis  and  0.1241/?^^^  N-m-s  for  the  reduced  dimensional 

state  space  analysis,  the  difference 
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is  expected. 


The 


1.0069/?^^^  N-m-s  sphere  is  based  only  on  neglecting  elliptic 

singularities  while  this  state  space  analysis  avoids  all 
singularities.  Reanalyzing  the  gimbal  limits  by  allowing 
hyperbolic  singularities  within  the  operating  momentum 
volume  will  increase  the  available  performance. 
Additionally,  the  gimbal  limits  are  based  on  finding  the 
largest  singularity-free  square  or  hypercube  due  to  the 
simplicity  of  implementing  the  search  algorithm  as  well  as 
the  simplicity  of  analyzing  the  resulting  gimbal  limits. 
This  provides  an  additional  performance  constraint  as  it 
does  not  include  all  contiguous  states  in  the  maximum  area 
or  hypervolume  that  are  singularity  free.  Finally,  a  voxel 
based  maximal  surface  detection  algorithm  may  improve  the 
performance  of  the  state  space  analysis  with  respect  to  the 
maximum  sphere  metric. 
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Gimbal  Limited  Momentum  Space  Comparison 
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Figure  45.  Four  CMG  singularity  free  momentum  space 

comparison 


Due  to  the  difficulty  in  visually  comparing  the 
complex  surfaces  of  singularities  [4,17]  and  the  complex 
shape  of  Figure  44,  further  study  is  required  to  verify 
that  the  results  of  the  state  space  analysis  matches  those 
found  in  the  literature. 

Additional  work  is  required  to  fully  determine 
whether  reduced  dimensional  state  space  analysis  provides 
sufficient  performance  or  if  a  full  dimensional  state  space 
analysis  is  required  to  maximize  array  performance. 
Ultimately,  the  result  may  depend  on  the  application. 
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Figure  46.  Maximum  sphere  superimposed  on  the  reduced 

dimensional  state  space  analysis-based  singularity 

free  momentum  space 
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Gimbal  Limited  Momentum  Space 
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Figure  47.  Maximum  sphere  superimposed  on  the  full 

dimensional  state  space  analysis-based  singularity 

free  momentum  space 


C .  REACTION  FORCES 

1 .  Validation  and  Verification 

Due  to  the  lack  of  material  in  literature  as  compared 
to  torque  and  momentum  space  studies,  rigorous  validation 
and  verification  of  the  reaction  force  equations  was 
required.  The  general  methodology  used  to  conduct  the 
validation  and  verification  involved  identifying  particular 
initial  conditions  that  output  a  simple,  intuitive  result, 
preferably  involving  a  single  parameter  such  as  axial 
torque  or  transverse  force. 
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There  are  two  primary  aspects  that  need  to  be 

verified.  The  first  is  the  inertial  motion  of  the  MED. 
This  is  the  right  hand  side  of  the  force-balance  equation. 
The  inertial  motion  equation  determines  the  linear  and 
rotational  trajectories  based  on  known  system  geometry  and 
simulated  or  real-world  telemetry.  The  other  aspect  is  the 
left  hand  side  of  the  force-balance  equation.  It  is  the 
summation  of  all  moments  and  forces  and  is  used  to 
determine  the  reaction  forces  and  moments.  Within  those 
two  primary  sections  are  a  total  of  six  parts  to  be 

verified:  forces  in  the  axial  and  both  transverse 

directions  and  moments  in  the  axial  and  both  transverse 
directions . 

a.  Reaction  Wheel  Inertial  Motion  Verification 

For  the  reaction  wheel  tests,  the  initial 

conditions  and  results  are  given  in  Table  4  through  Table 
9.  The  RW  is  assumed  to  be  oriented  such  that  the  spin 

axis  is  parallel  with  the  z-axis  of  the  spacecraft  and  has 
the  following  physical  characteristics: 

'1  0  O' 

0  1  0  kg-m^ 

0  0  1 

m  =  \  kg 

The  first  test  examines  rotations  in  the  axial 

direction  and  is  broken  up  into  four  parts.  The  first  two 

parts  have  the  RW  centered  in  the  spacecraft  body  and  have 

the  RW  accelerate  at  a  constant  rate  and  the  spacecraft 

accelerate  at  a  constant  rate,  respectively.  The  results 

for  these  tests  give  the  expected  inertial  moment  in  the 

axial  direction.  The  third  part  involves  the  RW 
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accelerating  at  a  constant  rate  in  one  direction  and  the 
spacecraft  accelerating  at  the  same  rate  in  the  opposite 
direction.  These  two  motions  cancel  each  other  out, 
resulting  in  no  inertial  moments  or  forces  applied  to  the 
RW.  The  final  part  has  the  RW  displaced  from  the  center  of 
the  spacecraft  and  has  the  spacecraft  rotating  at  a 
constant  rate.  The  result  shows  an  inertial  force  oriented 
towards  the  center  of  the  spacecraft. 


Applied  Ang.  Accel, 
(rad/sec^) 

Initial  Ang.  Rates  (rad/sec) 

RW  Location  (m) 

S/C 

RW 

S/C 

RW 

Xs/c 

Ys/c 

^^9 

X 

y 

z 

z 

mm 

n 

z 

z 

Testl 

a 

0 

0 

0 

0 

0 

0 

0.1 

0 

0 

0 

0 

b 

0 

0 

0 

0 

0 

0.1 

0 

0 

0 

0 

0 

c 

0 

0 

0 

0 

0 

-0.1 

0.1 

0 

0 

0 

0 

d 

0.1 

0 

0 

0 

0 

0 

0 

0 

0 

a 

0 

Table  4 .  RW  V&V  test  one  initial  conditions 


RW  Inertial  Force  (N) 

RW  Inertial  Moment  (N-m) 

X 

D 

z 

X 

y 

z 

a 

0 

0 

0 

0 

0 

0.1 

b 

0 

0 

0 

0 

0 

0.1 

c 

0 

0 

0 

0 

0 

0 

d 

-0.1 

0 

0 

0 

0 

0 

Table  5.  RW  V&V  test  one  results 


Applied  Ang.  Accel, 
(rad/sec^) 

Initial  Ang.  Rates  (rad/sec) 

RW  COM  Location  (m) 

S/C 

RW 

S/C 

RW 

X 

y 

z 

z 

mm 

n 

z 

z 

CM 

a 

0 

0 

0 

0 

0.1 

0 

0 

0 

0 

0 

0 

•M 

(/) 

0) 

b 

0 

0 

0 

0.1 

0 

0 

0 

0 

0 

0 

0 

H 

c 

0 

0 

0 

0 

0.1 

0 

0 

0 

0 

0 

0.2 

Table  6.  RW  V&V  test  two  initial  conditions 
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RW  Inertial  Force  (N) 

RW  Inertial  Moment  (N-m) 

X 

y 

z 

X 

Y 

z 

CM 

a 

0 

0 

0 

0 

0.1 

0 

•M 

01 

V 

0 

0 

0 

0.1 

0 

0 

H 

c 

0 

0 

0 

* 

0.1 

0 

Table  7 .  RW  V&V  test  two  results 


* 


5 


-m) 


The  second  test  is  composed  of  three  parts  and 
examines  transverse  motions.  For  all  parts  of  this  test, 
the  RW  is  located  at  the  center  of  the  spacecraft.  Parts 
one  and  two  have  the  spacecraft  with  a  rotational 
acceleration  along  a  single  transverse  axis.  The  results 
show  an  inertial  moment  along  the  corresponding  axis  with 
the  appropriate  magnitude.  The  third  part  repeats  the 
first,  however  the  RW  is  now  rotating  at  a  constant  rate 
along  its  spin  axis.  The  results  show  a  constant  inertial 
moment  along  the  appropriate  transverse  axis  as  well  as  a 
linearly  increasing  inertial  moment  about  the  other 
transverse  axis.  This  increasing  moment  matches  the 
expected  behavior  of  a  CMG,  which  is  what  the  initial 
conditions  for  the  third  part  simulated. 


Applied  Linear  Accel,  (m/sec^) 

RW  COM  Location  (m) 

S/C 

Xs/c 

Ys/c 

Zs/c 

X 

y 

z 

no 

a 

0 

0 

0 

0.1 

0 

0 

•M 

iA 

Q) 

b 

0 

0 

0 

0 

0.1 

0 

H 

c 

0 

0 

0 

0 

0 

0.1 

Table  8 .  RW  V&V  test  three  initial  conditions 
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RW  Inertial  Force  (N) 

RW  Inertial  Moment  (N-m) 

X 

y 

z 

X 

y 

z 

m 

a 

0.1 

0 

0 

0 

0 

0 

•M 

iA 

01 

V 

0 

0.1 

0 

0 

0 

0 

H 

c 

0 

0 

0.1 

0 

0 

0 

Table  9.  RW  V&V  test  three  results 


The  final  test  involved  linear  acceleration  along 
all  three  RW  axes.  The  results  all  showed  an  inertial 
force  with  the  appropriate  magnitude  in  the  appropriate 
direction . 


b.  RW  Transverse  Reaction  Verification 

The  reaction  verification  is  broken  up  into  two 
parts:  analyzing  an  applied  inertial  transverse  force  and 
an  applied  inertial  transverse  moment.  Both  the  moment  and 
force  are  applied  to  act  on  the  center  of  mass  (COM) .  Each 
part  is  then  subdivided  into  four  sections.  The  first 
section  analyzes  the  effect  of  moving  the  COM  from  one 
bearing  to  the  other,  along  the  length  of  the  shaft.  The 
second  section  analyzes  the  effect  of  shaft  length.  Both 
sections  are  performed  for  both  the  cantilever  and  fixed- 
fixed  end  condition  cases.  The  RW  is  assumed  to  have  the 
same  physical  properties  as  described  in  the  inertial 
motion  verification.  Additionally,  the  nominal  shaft 

length  is  one  meter  and  the  nominal  COM  location  is  at  the 
midpoint  of  the  shaft.  The  applied  force  has  a  magnitude 
of  1  N  and  the  applied  moment  has  a  magnitude  of  1  N-m. 
These  values  were  chosen  to  generate  plots  whose 
relationships  could  be  scaled  to  fit  the  scale  of  the  MED 
array  being  analyzed. 
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Figure  49  through  Figure  52  give  the  transverse 
force  case  with  Figure  48  giving  the  direction  of  the 
applied  force.  In  accordance  with  the  convention  described 
in  Chapter  V  Section  B,  bearing  one  is  at  the  bottom  of  the 
RW  and  bearing  two  is  at  the  top. 
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Figure  48.  Applied  inertial  force  geometry  for 

RW  reaction  validation 
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COM  Location  Variation 


Figure  49.  Fixed-fixed  RW  COM  variation  reaction 
validation  for  an  applied  inertial  force 
at  the  center  of  the  shaft 
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Figure  50.  Fixed-fixed  RW  shaft  length  variation 

reaction  validation  for  an  applied  inertial 
force  at  the  center  of  the  shaft 
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COM  Location  Variation 
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Figure  51.  Cantilever  RW  COM  variation  reaction 
validation  for  an  applied  inertial  force 
at  the  center  of  the  shaft 


RW  Shaft  Length  Variation 


Figure  52.  Cantilever  RW  shaft  length  variation 
reaction  validation  for  an  applied 
inertial  force  at  the  center  of  the  shaft 
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For  both  the  fixed-fixed  and  cantilever  cases, 
the  reaction  forces  and  moments  behave  as  expected  with 
respect  to  the  geometrical  variations. 

Figure  53  shows  the  geometry  for  the  applied 
moment  section,  with  Figure  54  to  Figure  57  giving  the 
results . 
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Figure  53.  Applied  inertial  moment  geometry  for 

RW  reaction  validation 
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COM  Location  Variation 


Figure  54.  Fixed-fixed  RW  COM  variation  reaction 
validation  for  an  applied  inertial  moment 
at  the  midpoint  of  the  shaft 
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Figure  55.  Fixed-fixed  RW  shaft  length  variation 

reaction  validation  for  an  applied  inertial 
moment  at  the  midpoint  of  the  shaft 
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COM  Location  Variation 
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Figure  56.  Cantilever  RW  COM  variation  reaction 
validation  for  an  applied  inertial  moment 
at  the  midpoint  of  the  shaft 
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Figure  57.  Cantilever  RW  shaft  length  variation 

reaction  validation  for  an  applied  inertial 
moment  at  the  midpoint  of  the  shaft 
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All  reaction  forces  and  moments  acted  as  expected 
through  this  analysis,  which  validate  the  application  of 
the  beam  tables. 

c.  Reaction  Wheel  Axial  Reaction  Verification 

Due  to  the  assumed  construction  of  the  fixed- 
fixed  system,  bearing  one  supports  all  axial  force  and 
moment  loads.  This  matches  the  reaction  of  the  cantilever 
system.  Since  Equations  (71)  and  (81)  explicitly  enforce 
this  property,  additional  analysis  is  not  required. 

d.  CMG  Inertial  Motion  Verification 

As  the  inertial  motion  of  the  gimbal  frame 

follows  the  same  kinematics  as  the  RW  frame,  the  prior 
analysis  sufficiently  proves  the  equations  used.  The  rotor 
frame,  however,  has  the  added  complexity  of  gimbal  motion 
and  so  must  be  analyzed. 

Due  to  the  complexity  of  the  resulting  inertial 

forces  and  moments  on  the  rotor  from  combined  gimbal  and 
spacecraft  rotations,  the  utility  of  performing  as 
exhaustive  an  analysis  as  was  performed  on  the  RW  equations 
is  limited.  Therefore,  two  specific  rotor  cases  will  be 

examined:  motion  by  the  gimbal  only  and  motion  by  the 
gimbal  with  counter-motion  by  the  spacecraft  (similar  to 
test  l.c  for  the  RW)  .  The  mechanical  properties  of  the 

rotor  will  match  those  of  the  RW,  the  origin  of  the  gimbal 
frame  will  be  coincident  with  the  origin  of  the  spacecraft 
body  frame,  and  the  gimbal  rotation  axis  will  be  parallel 
to  the  Z-axis  of  the  spacecraft  body.  All  bodies  will 
start  at  rest.  Table  10  gives  the  specific  conditions  of 
the  test. 
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Applied  Ang.  Accel,  (rad/sec^) 

Rotor  COM  Location  (m) 

S/C 

Gimbal 

Xg 

Yg 

Zg 

X 

y 

z 

S 

Testl 

0 

0 

0 

0 

0 

0 

0.1 

Test  2 

0 

0 

0 

0 

0 

-0.1 

0.1 

Table  10.  CMG  V&V  initial  conditions 


Figure  58  and  Figure  59  give  the  results  of  the 
tests.  The  small  constant  positive  inertial  force  in  the  y 
direction  in  test  one  is  caused  by  the  angular  acceleration 
of  the  gimbal .  The  decreasing  force  in  the  negative  x 
direction  represents  the  centrifugal  force  caused  by  the 
actual  rotation  of  the  gimbal.  As  the  gimbal  rate 
increases,  the  magnitude  of  the  force  increases  as  seen  in 
Equation  (106)  . 

" =  V,  +  :R 

+2 IR (10  6) 
"a,  =  + 
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Figure  58.  Results  of  the  first  CMG  V&V  test 


io-’5  Inertial  Motion  CMG  Test  2 
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Figure  59.  Results  of  the  second  CMG  V&V  test 
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For  the  second  test,  the  prominent  feature  is  the 
noise  in  the  inertial  force  plot.  Though  of  a  negligible 
magnitude,  it  likely  appears  due  to  several  adverse  factors 
including  the  rudimentary  simulation  method  and  relatively 
rapidly  rotating  bodies.  The  important  feature  is  the  fact 
that  the  counter-spinning  bodies  remove  the  inertial  moment 
from  the  rotor  as  it  now  appears  to  stand  still  in  the 
inertial  frame. 

2 .  Reaction  Wheel  Maneuvers 

The  spacecraft  geometry  for  the  RW  analysis  is  given 
in  Figure  60  and  tabulated  in  Table  11.  The  maneuver  is  an 
arbitrary  rest  to  rest  slew  from  one  orientation  to 
another.  The  analyzed  telemetry  data  includes  only  the 
spacecraft  quaternions  and  body  rates,  the  RW  angular 
acceleration  and  an  array  of  the  time  steps. 
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Figure  60.  Geometry  detail  for  the  RW  maneuver  force 

and  moment  analysis 
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RW 

r  (deg) 

J3  (deg) 

COM  (Xb/Yb/Zb)  (m) 

1 

0 

54.73 

0.25 

0 

0 

2 

90 

54.73 

0 

0.25 

0 

3 

180 

54.73 

-0.25 

0 

0 

4 

270 

54.73 

0 

-0.25 

0 

Table  11.  Tabulated  values  for  RW  array  layout  for 

maneuver  analysis 

The  RWs  used  in  the  analysis  have  the  same  properties 
as  those  used  in  the  momentum  space  case  study.  Telemetry 
was  derived  from  the  model  developed  in  [4] . 

a.  Eigenaxis  Maneuver 

The  first  maneuver  is  based  on  an  eigenaxis 
controller  where  the  eigenaxis  between  the  initial  and 
final  state  is  calculated  and  a  torque  command  is  developed 
to  rotate  the  spacecraft  about  this  axis.  Figure  61 
describes  the  basic  telemetry  for  the  maneuver. 


ill 


Wheel  Rates 


Body  Rate  Norm 


Figure  61.  Eigenaxis  slew  maneuver  for  RW  array  force 

and  moment  analysis 

During  the  maneuver,  the  first  and  third  RW 
saturate  at  +/-  4000  RPM  to  achieve  a  constant  spacecraft 
body  rate  magnitude  of  slightly  above  4  deg/sec  about  the 
eigenaxis.  The  calculated  inertial  force  and  moment 
magnitudes  are  given  in  Figure  62  and  Figure  63. 
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Inertial  Force  Magnitude 


Figure  62.  Inertial  force  magnitudes  during  an  eigenaxis 

slew  maneuver 


Inertial  Moment  Magnitude 
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Figure  63.  Inertial  moment  magnitudes  during  an  eigenaxis 

slew  maneuver 


During  this  maneuver,  the  magnitudes 
inertial  forces  are  equal  for  opposite  pairs  of 


of  the 
(one 
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RWs 


and  three,  two  and  four) .  This  matches  expectations  as  the 
RWs  lie  equidistant  along  a  common  axis  from  the  spacecraft 
body  origin  which  is  the  point  about  which  the  spacecraft 
rotates.  Thus,  one  would  expect  the  centrifugal  force  on 
the  RWs  to  be  equal  in  magnitude  and  opposite  in  direction. 
The  difference  in  magnitude  between  the  two  pairs  indicates 
that  RWs  two  and  four  are  moving  faster  through  inertial 
space,  indicating  that  the  projection  of  the  eigenaxis 
along  the  line  connecting  RW  two  and  RW  four  is  larger  than 
the  projection  on  the  line  connecting  the  other  set.  This 
assumption  is  borne  out  by  the  magnitude  of  the  inertial 
moments.  Since  the  spacecraft  is  spinning  about  an  axis 
that  is  closer  to  the  rotation  axis  of  RW  one  and  three, 
they  will  have  a  larger  inertial  moment  induced  on  them. 

The  steps  seen  at  about  12  seconds  and  46  seconds 
indicate  the  times  when  the  RW  angular  acceleration  goes 
from  non-zero  to  zero  (and  vice-versa) . 

The  reaction  force  plots  for  the  fixed-fixed 
configuration  given  in  Figure  64  and  Figure  65  show  an 
interesting  property.  The  inertial  moments  are  between  one 
and  two  orders  of  magnitude  larger  than  the  forces.  From 
Figure  54,  a  one  meter  shaft  results  in  reaction  forces 
about  1.5  times  larger  than  the  applied  moment. 
Furthermore,  Figure  55  shows  that  a  0.2  m  shaft,  the 
assumption  used  for  this  analysis,  results  in  reaction 
forces  eight  times  larger  than  the  applied  moment.  The 
inertial  moment,  therefore,  contributes  most  to  the 
reaction  force  at  the  bearings. 
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BRG1  Rxn  Forces 


Figure  64 .  Bearing  one  reaction  force  magnitudes  for  the 
eigenaxis  slew,  fixed-fixed  configuration 


BRG2  Rxn  Forces 


Figure  65.  Bearing  two  reaction  force  magnitudes  for  the 
eigenaxis  slew,  fixed-fixed  configuration 


That  fact  explains  why  both  bearings  have 
approximately  the  same  force  magnitude,  despite  the  fact 


115 


that  bearing  two  does  not  support  a  force  in  the  axial 
direction.  When  examining  the  reaction  moments.  Figure  66 
and  Figure  67,  bearing  one  shows  a  significant  larger 
moment  than  bearing  two  during  wheel  operation.  This  is 
attributed  to  the  fact  that  only  bearing  one  resists  an 
axial  moment  due  to  the  motor  being  co-located  with  bearing 
one.  The  larger  reaction  moment  magnitude  for  RWs  one  and 
three  stems  from  the  fact  that  they  are  accelerating  at  a 
higher  rate. 

When  the  RWs  are  not  accelerating,  the  reaction 
moment  magnitudes  on  both  bearings  roughly  match  each 
other . 


BRG1  Rxn  Moments 


Figure  66.  Bearing  one  reaction  moment  magnitudes  for  the 
eigenaxis  slew,  fixed-fixed  configuration 
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BRG2  Rxn  Moments 


Figure  67 .  Bearing  two  reaction  moment  magnitudes  for  the 
eigenaxis  slew,  fixed-fixed  configuration 

With  the  cantilever  case,  the  bearing  one 
reaction  forces  and  moments  act  as  expected,  shown  in 
Figure  68  and  Figure  69.  The  small  inertial  force 
magnitude  and  short  lever  arm  prevent  the  creation  of  a 
noticeable  reaction  moment  at  the  bearing,  thus  the 
reaction  moment  and  inertial  moment  plots  appear  to  match. 

The  magnitude  of  the  reaction  force  on  bearing 
one  in  the  cantilever  case  is  much  less  than  for  the  fixed- 
fixed  case  because,  according  to  the  standard  beam  tables 
[13,  pp .  208-223],  the  inertial  moment  does  not  result  in  a 
reaction  force  on  the  bearings. 
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BRG1  Rxn  Forces 


Figure  68.  Bearing  one  reaction  force  magnitudes  for  the 
eigenaxis  slew,  cantilever  configuration 
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Figure  69.  Bearing  one  reaction  moment  magnitudes  for  the 
eigenaxis  slew,  cantilever  configuration 
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b.  Optimal  Maneuver 


The  second  maneuver  analyzed  is  an  optimal 
maneuver  where  the  cost  function,  J  ,  represents  a  quantity 
that  the  controller  minimizes.  This  cost  function  can 
include  such  quantities  as  elapsed  time  or  wheel  power. 
The  same  RW  configuration  is  used  and  the  telemetry  is 
derived  from  the  same  source.  The  maneuver  is  conducted 
with  the  same  initial  and  final  conditions  as  before. 
Figure  70  shows  the  maneuver. 


Wheel  Rates 


Figure  70.  Optimal  slew  maneuver  for  RW  array  force  and 

moment  analysis 

For  this  maneuver,  the  wheel  rates  of  all  wheels  are 
now  closer  to  or  at  saturation  for  a  longer  period  of  time. 
The  body  rate  magnitude  is  somewhat  higher  throughout  the 
maneuver,  and  the  final  state  is  reached  slightly  sooner 
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than  the  eigenaxis  maneuver.  This  additional  exercising  of 
the  RWs  is  typical  of  an  optimal  maneuver  as  the  algorithm 
seeks  to  use  the  excess  performance  that  is  present  but 
that  the  eigenaxis  slew  ignores. 

Since  there  is  no  single  axis  of  rotation  for  the 
maneuver  and  the  RWs  tend  to  continuously  accelerate  or 
decelerate,  the  inertial  force  and  moment  magnitudes 
induced  on  the  RWs  tend  to  be  more  complicated,  as  shown  in 
Figure  71  and  Figure  72. 


Inertial  Force  Magnitude 


Figure  71.  Inertial  force  magnitudes  during  an  optimal 

slew  maneuver 
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Inertial  Moment  Magnitude 


Figure  72.  Inertial  moment  magnitudes  during  an  optimal 

slew  maneuver 

There  are  some  similarities  between  the  time  optimal 
and  eigenaxis  slew  maneuvers  with  regards  to  the  inertial 
forces  and  moments.  As  with  the  eigenaxis  slew,  the 
optimal  slew  produces  inertial  forces  that  are  equal  in 
magnitude  to  the  same  paired  RWs :  one  and  three,  two  and 
four.  This  equality  is  expected  due  to  the  positional 
symmetry  in  the  system.  Since  all  RWs  are  equidistant  from 
the  z-axis,  a  spacecraft  body  rotation  about  this  axis 
produces  equal  inertial  force  in  all  components.  RW  one 
and  three  are  positioned  along  the  x-axis  equidistant  from 
the  origin,  so  a  spacecraft  body  rotation  rate  about  the  y- 
axis  produces  equal  and  opposite  inertial  forces  in  that 
pair,  and  zero  inertial  force  in  the  other.  The  opposite 
is  true  for  a  spacecraft  body  rotation  rate  about  the  x- 
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axis.  The  inertial  forces  tend  to  be  slightly  larger  for 
all  four  components  due  to  the  higher  body  rates;  however 
they  are  still  an  order  of  magnitude  below  the  inertial 
moments . 

When  analyzing  the  behavior  of  the  moments,  though,  it 
is  possible  to  determine  when  changes  in  wheel  acceleration 
occur  as  these  accelerations  cause  moments  that  tend  to 
dominate.  This  is  expected  as  the  wheel  can  change 
acceleration  much  faster  than  the  spacecraft  will  alter  its 
rotational  motion.  Removing  the  torque  induced  by  the 
wheel  itself  would  produce  much  smoother  plots,  as  seen 
during  the  periods  of  constant  wheel  acceleration  in  Figure 
63. 

The  bearing  force  and  moment  reaction  magnitudes  are 
given  in  Figure  73  through  Figure  76.  As  before,  the 
inertial  moments  dominate  the  reaction  forces,  and  bearing 
one  supports  much  larger  reaction  moments  than  bearing  two 
due  to  the  actual  torque  produced  by  wheel  accelerations. 
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BRG1  Rxn  Forces 


Figure  73.  Bearing  one  reaction  force  magnitudes  for  an 
optimal  slew,  fixed-fixed  configuration 


BRG2  Rxn  Forces 


Figure  74.  Bearing  two  reaction  force  magnitudes  for  an 
optimal  slew,  fixed-fixed  configuration 
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BRG1  Rxn  Moments 


Figure  75.  Bearing  one  reaction  moment  magnitudes  for  an 
optimal  slew,  fixed-fixed  configuration 


BRG2  Rxn  Moments 


Figure  76.  Bearing  two  reaction  moment  magnitudes 

for  an  optimal  slew,  fixed-fixed  configuration 
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Figure  77  and  Figure  78  give  the  cantilever  results. 
The  reaction  forces  and  moments  match  the  inertial  forces 
and  moments,  an  expected  result. 

Although  the  optimal  maneuver  results  in  more  active 
RW  utilization,  the  force  and  moment  magnitudes  seen  by  the 
bearings  are  of  roughly  equivalent  magnitudes  to  those  seen 
during  the  eigenaxis  maneuver.  Because  the  wheel 
accelerations  tend  to  dominate  the  inertial  moment  equation 
(and  the  inertial  moment  dominates  the  inertial  force  with 
respect  to  bearing  reactions) ,  the  optimal  maneuver  tends 
to  maintain  a  higher  average  loading.  This  particular 
eigenaxis  maneuver,  however,  produces  very  large  and  rapid 
changes  in  force  and  moment  magnitudes  due  to  the  bang-off- 
bang  nature  of  the  maneuver.  This  shock  to  the  system  is 
not  present  in  the  optimal  maneuver. 


BRG1  Rxn  Forces 


Figure  77.  Bearing  one  reaction  force  magnitudes  for  an 
optimal  slew,  cantilever  configuration 
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BRG1  Rxn  Moments 


Figure  78.  Bearing  one  reaction  moment  magnitudes  for  an 
optimal  slew,  cantilever  configuration 

3 .  CMG  Maneuvers 

Using  four  CMGs  arranged  in  the  same  geometrical  shape 
shown  in  Figure  60  and  described  in  Table  11,  reaction 
analyses  for  an  eigenaxis  slew  and  an  optimal  slew  were 
performed.  In  this  case,  the  initial  and  final  conditions 
between  the  two  maneuvers  are  slightly  different  because 
the  telemetry  data  is  from  a  higher  fidelity  simulation  of 
an  imaging  spacecraft.  The  CMG  mechanical  properties  are 
the  same  as  those  for  the  momentum  and  torque  space 
analysis  in  Chapter  VI  Sections  A  and  B  with  the  addition 
that  the  rotor  shaft  length  is  assumed  to  be  0.1  m,  the 
gimbal  shaft  length  is  assumed  to  be  0.5  m  and  bearing  one 
for  the  rotor  is  assumed  to  be  mounted  0.05  m  from  the 
gimbal  COM  in  the  rotor's  axial  direction  (x^=x^).  The  CMG 
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has  an  assumed  mass  of  2  kg.  The  assumed  rotor  inertia 
tensor  is: 

'.0429  0  0 

7^^=  0  .0429  0  kg-rn^  . 

0  0  .0025 

Although  there  are  four  possible  variations  for  the 
gimbal  and  rotor  end  point  conditions,  only  the  fixed- 
f ixed/f ixed-f ixed  case  will  be  examined  as  it  provides 
sufficient  detail  for  this  analysis.  A  more  rigorous 

examination  would  be  warranted  in  other  cases  such  as 
actual  satellite  or  CMG  mechanism  design. 

a.  Eigenaxis  Maneuver 

The  analyzed  maneuver  is  actually  the  end  of  one 
maneuver  and  the  beginning  of  another.  Due  to  the  nature 
of  the  data,  representing  an  entire  maneuver  is  impractical 
due  to  the  length  of  time  involved.  Analyzing  the  impulses 
between  two  maneuvers  reduces  the  amount  of  time  to  be 
analyzed  while  still  examining  two  major  control  events. 

Initially  the  spacecraft  is  rotating  at  a  small 
rate.  The  CMGs  are  gimbaled  to  reduce  the  body  rate  for  a 
few  seconds .  They  are  then  gimbaled  again  to  induce  a  new 
spacecraft  body  rate.  Figure  79  through  Figure  85  give  the 
maneuver  telemetry  and  results  of  the  force  analysis.  The 
large  force  of  36  N  seen  at  41  and  45  seconds  in  the  first 
gimbal  bearing  is  caused  by  the  inertial  moment  induced  on 
the  rotor  at  those  times.  Due  to  the  complexity  of  the 
system  dynamics,  this  result  cannot  be  easily  gleaned  from 
the  analysis  of  the  equations  of  motion  as  was  done  in 
Equation  (106)  . 
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Due  to  the  short  length  of  the  rotor  shaft,  the 
reaction  forces  are  many  times  larger  than  the  inertial 
moment  would  imply.  The  reaction  force  at  the  rotor 
bearing  one  is  then  transmitted  to  the  gimbal  in  the  axial 
direction.  Since  the  second  gimbal  bearing  does  not 

support  forces  in  this  direction,  the  entire  load  is  taken 
up  by  the  first  bearing  (or  more  specifically,  the 
gimbaling  motor  located  at  that  bearing) . 


Body  Rate  Norm 


Figure  79.  Eigenaxis  slew  telemetry  for  CMC  force  analysis 
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^-,0-3  Inertial  Force  Magnitude 


Figure  80.  Inertial  forces  from  the  CMG  force  analysis 

eigenaxis  slew 


Inertial  Moment  Magnitude 


Figure  81.  Inertial  moments  from  the  CMG  force  analysis 

eigenaxis  slew 
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BRG1  Rxn  Forces 


Figure  82  . 


Bearing  one  reaction  force  magnitudes  from  the 

eigenaxis  slew 


BRG2  Rxn  Forces 


Figure  83.  Bearing  two  reaction  force  magnitudes  from  the 

eigenaxis  slew 
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BRG1  Rxn  Moments 


Figure  84.  Bearing  one  reaction  moment  magnitudes  from  the 

eigenaxis  slew 


BRG2  Rxn  Moments 


Figure  85.  Bearing  two  reaction  moment  magnitudes  from  the 

eigenaxis  slew 
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Although  the  maneuvers  between  the  RW  and  CMG 
analyses  are  very  different,  it  appears  that  a  CMG  has  the 
potential  to  induce  much  larger  forces  on  its  bearings 
compared  to  a  RW.  This  is  also  how  the  CMG  provides  large 
torques  to  the  spacecraft.  Thus,  the  bearings  for  a  CMG 
must  be  built  to  withstand  larger  loads  than  those  for  a 
comparably  sized  RW. 

b.  Optimal  Maneuver 

Similar  to  the  eigenaxis  maneuver,  the  optimal 
maneuver  shows  the  completion  and  commencement  of  two  back- 
to-back  slews.  The  gimbal  accelerations  are  not  as 
excessive  in  this  case  which  slightly  reduces  the  induced 
rotor  inertial  moment.  This,  in  turn,  reduces  the  large 
reaction  force  seen  at  the  gimbal/rotor  junction  and  the 
first  gimbal  bearing.  The  results  are  given  in  Figure  86 
through  Figure  92 . 
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Gimbal  Rates 


o 


Figure  86.  Optimal  slew  telemetry  for  CMG  force  analysis 


Inertial  Force  Magnitude 


Figure  87.  Inertial  force  magnitudes  from  the  CMG  force 

analysis  optimal  slew 
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Figure 


Figure 


^0-3  Inertial  Moment  Magnitude 


88.  Inertial  moment  magnitudes  from  the  CMG  force 

analysis  optimal  slew 


BRG1  Rxn  Forces 


89.  Bearing  one  reaction  force  magnitudes  for  the 

optimal  slew 
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BRG2  Rxn  Forces 


Figure  90.  Bearing  two  reaction  force  magnitudes  for  the 

optimal  slew 


BRG1  Rxn  Moments 


Figure  91.  Bearing  one  reaction  moment  magnitudes  for  the 

optimal  slew 
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BRG2  Rxn  Moments 


Figure  92 .  Bearing  two  reaction  moment  magnitudes  for  the 

optimal  slew 

Contrary  to  what  was  seen  for  the  RWs  where  the 
peak  measurements  between  the  control  cases  were  roughly 
equal,  the  peak  reaction  forces  for  the  eigenaxis  maneuvers 
are  slightly  larger  than  those  seen  during  the  optimal 
maneuvers.  The  optimal  maneuver,  however,  maintains  the 
reaction  force  for  a  larger  amount  of  time  while  the 
eigenaxis  maneuver  introduces  more  abrupt  force  changes. 
The  same  trend  is  seen  for  the  reaction  moments.  This 
indicates  that,  for  this  particular  maneuver  with  this 
particular  geometry,  the  optimal  controller  is  easier  on 
the  components  than  the  eigenaxis  maneuver. 
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D.  SUMMARY 

Chapter  VI  presented  a  series  of  applications  of  the 
equations  and  algorithms  derived  in  the  prior  chapters. 
Overall,  the  algorithms  were  shown  to  be  able  to  rapidly 
analyze  many  different  variations  on  a  potential  ACS,  as 
well  as  the  ability  to  analyze  maneuvers  to  determine 
expected  reaction  forces  and  moments  using  only  MED 
geometry  and  standard  telemetry  data,  allowing  for 
comparison  of  different  controllers  and  different  array 
geometries . 

1 .  Momentvim  Envelope 

For  the  RW  momentum  envelopes,  the  equations  were 
validated  by  analyzing  the  geometry  of  the  momentum  spaces 
with  the  geometry  of  different  arrays.  Additionally,  an 
analysis  was  performed  on  the  effect  of  rotor  imbalance  and 
indicated  that,  for  expected  imbalance  magnitudes,  there 
was  essentially  no  performance  loss  in  the  array. 

The  CMC  MIA  analysis  was  effective  in  rapidly 
generating  a  highly  detailed  maximum  momentum  space  for 
different  CMC  array  geometries.  The  associated  angle 
residual  plots  gave  an  indication  where  interior 
singularity  surfaces  met  the  maximum  surface. 

Finally,  the  algorithms  were  used  to  generate  a 
maximum  momentum  envelope  for  a  hybrid  RW/CMG  array  by 
summing  the  individually  generated  RW  array  momentum 
envelopes  and  CMC  array  envelopes.  This  resulted  in  the 
possibility  of  estimating  the  hybrid  array  performance  with 
respect  to  spacecraft  slews  as  well  as  introducing  the  idea 
of  adding  RWs  to  a  CMC  array  to  allow  the  CMGs  to  steer 
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allowing  expanded 


around  elliptic  singularities  thus 
utilization  of  the  CMG  performance  envelope. 

2 .  Torque  Envelopes 

The  torque  envelope  was  shown  to  be  of  the  same  shape 
as  the  momentum  envelope  for  RW  arrays. 

For  CMG  arrays,  the  state  space  analysis  algorithms 
were  implemented  and  used  to  generate  singularity  free 
momentum  spaces.  It  was  discovered  that  the  current  state 
space  analysis  methods  resulted  in  a  momentum  envelope  that 
allowed  about  12%  of  the  capability  of  current  CMG  array 
capability  for  the  semi-adjacent  analysis,  and  about  6%  for 
the  fully  adjacent  analysis.  This  had  several  causes,  the 
primary  one  being  that  the  current  performance  metric  is 
based  on  including  hyperbolic  singularities  while  the  state 
space  analyses  excluded  all  singularities.  A  theoretical 
control  algorithm  was  presented,  and  an  initial  test  was 
performed  showing  that  both  the  control  algorithm  concept 
was  valid  and  the  momentum  space  held  no  singularities 
along  the  trajectory. 

3.  Reaction  Force  Analysis 

The  developed  equations  were  validated  for  RWs  and 
CMGs  using  simple  test  cases.  Third  party  simulators  were 
used  to  generate  telemetry  to  analyze  different 
controllers,  analyze  different  RW  and  CMG  construction 
assumptions  and  prove  the  feasibility  of  using  third-party 
telemetry  for  the  analysis. 

For  both  RW  and  CMG  arrays,  an  eigenaxis  slew  and  an 
optimal  slew  were  analyzed.  For  the  RW  analyses,  peak 

force  and  moment  magnitudes  were  roughly  equal.  The 
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optimal  slew  operated  the  RWs  at  a  higher  duty  cycle,  but 
provided  a  much  smoother  reaction  force  and  moment  history. 
The  eigenaxis  slew  tended  to  create  abrupt  changes  in  the 
reaction  forces  seen  on  the  system. 

For  the  CMG  arrays,  the  peak  magnitudes  for  the 
optimal  maneuver  were  less  than  for  the  eigenaxis  maneuver. 
The  smoothness  seen  in  the  RW  analysis  carried  over  to  the 
CMGs  for  the  optimal  maneuver,  and  the  abrupt  applications 
of  relatively  large  forces  and  moments  was  even  larger  for 
the  CMG  case  than  the  RW  case. 
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VII .  CONCLUSIONS  AND  FUTURE  WORK 


A.  CONCLUSIONS 

The  goal  of  this  thesis  was  to  develop  the  equations 
and  algorithms  necessary  to  automate  the  analysis  of  MED 
arrays  for  the  purpose  of  ACS  design  and  analysis.  This 
involved  looking  at  three  specific  areas:  momentum  envelope 
analysis,  torque  envelope  analysis,  and  reaction  force  and 
moment  analysis. 

For  RWs,  a  simple  matrix  equation.  Equation  (23) ,  was 
developed  that  can  rapidly  create  a  point  cloud  of  momentum 
states.  This  provided  a  means  to  estimate  possible  ACS 
performance  based  on  the  minimum  momentum  point  on  the 
maximum  surface.  This,  in  turn,  allowed  array  geometry  or 
the  physical  characteristics  of  the  RW  (such  as  inertia 
tensor  or  maximum  wheel  speed)  to  be  varied  and  the  result 
of  this  variance  with  respect  to  array  performance  to  be 
analyzed.  It  was  also  shown  that  the  RW  torque  space  has 

the  same  shape  as  the  RW  momentum  space,  with  the  torque 

space  magnitude  being  defined  by  the  maximum  wheel 
accelerations . 

For  CMGs,  the  MIA  process  was  developed  to  rapidly 
analyze  the  maximum  momentum  surface  by  removing  the  need 
to  perform  a  sweep  analysis  through  all  gimbal  angle 
combinations.  It  also  provided  an  indication  of  where 

interior  singular  surfaces  meet  the  maximal  momentum 

surface  via  an  angle  residual  plot. 

When  analyzing  torque  space,  it  was  shown  that  CMC 

arrays  produce  torque  that  is  state  dependent.  This  led  to 

the  development  of  the  state  space  analysis  in  order  to 
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generate  singularity  free  momentum  spaces  by  limiting  the 
gimbal  angles  of  each  CMG  in  the  array.  The  idea  of 
contiguousness  was  introduced  to  ensure  that  an  enclosed 
area  within  the  reduced  dimensional  state  space  analysis 
map  would  only  contain  states  that  yield  a  realizable 
gimbal  angle  trajectory.  Both  of  the  two  primary  state 
space  analyses,  reduced  dimensional  and  full  dimensional, 
successfully  created  singularity  free  spaces. 

The  reaction  force  and  moment  analysis  involved  a 
topic  that  had  received  little  discussion  in  the  researched 
literature.  The  equations  were  developed  using  a 
combination  of  first  principle  derivations  and  standard 
reference  tables.  The  analysis  was  formulated  such  that 
common  telemetry  data  and  known  MED  array  geometry, 
consistent  with  the  previous  momentum  and  torque  analyses, 
would  be  sufficient  to  solve  the  equations.  Results  from 
the  verification  and  validation  tests  showed  expected 
trends  and  provided  useful  relationships  between  MED 
geometry  changes  and  resulting  reaction  force/moment 
changes.  Initial  telemetry  testing  indicated  that  optimal 
controllers  and  eigenaxis  controllers  created  reaction 
forces  and  moments  that  were  of  the  same  magnitude  for  RWs, 
while  for  CMGs  the  optimal  controller  showed  smaller  force 
and  moment  magnitudes.  Additionally,  eigenaxis  controllers 
tended  to  impact  the  system  with  sudden,  large  forces  and 
moments  while  the  optimal  controllers  resulted  in  smoother 
reaction  forces  and  moments.  In  all  cases,  the  average 
duty  cycle  of  all  MEDs  operating  with  the  optimal 
controllers  was  higher  than  the  average  duty  cycle  under 
the  eigenaxis  controllers. 
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Overall,  this  thesis  showed  that  automated  analyses 
for  generic  MED  arrays  was  feasible,  and  should  prove 
useful  for  the  rapid  iteration  of  initial  ACS  design. 


B .  FUTURE  WORK 


Although  the  goal  of  developing  algorithms  to  rapidly 
analyze  arrays  of  RWs  and  CMGs  was  achieved,  there  are 
several  additional  features  that  can  be  added  to  improve 
analysis  fidelity  and  usefulness. 


1 .  Momentvim  Analysis 

•  Very  tight  tolerances  in  fly-wheel  balancing  are 
required  to  ensure  the  long  life  of  RWs  and  CMGs. 
This  does  allow  for  some  imbalances,  however,  and 
could  impact  the  performance  of  platforms  that 
require  very  tight  pointing  requirements. 
Especially  for  RWs,  it  would  be  useful  to  modify 
Equation  (28)  so  that  it  included  tracking  of  RW 
frame  rotation  as  the  wheel  rotates. 

•  For  CMGs,  the  rotor  provides  the  bulk  of  the 
component's  momentum.  The  rotation  of  the  gimbal, 
however,  adds  some  momentum  and  this  should  be 
included  as  part  of  the  analysis.  This  may  also 
benefit  the  singularity  free  state  space  analysis 
as  it  would  then  provide  for  the  possibility  of 
angular  momentum  vectors  outside  the  plane  of 
rotation . 

•  The  maximal  surface  detection  algorithm  is  a  very 
rudimentary  one  that  sections  a  unit  sphere  and 
finds  the  point  furthest  from  the  point  in  any 
given  section.  Anomalies  can  present  themselves 
at  the  north  and  south  poles  of  the  sphere.  A 
voxel-based  approach  could  provide  a  more  robust 
solution  as  well  as  a  possible  performance 
improvement . 

•  Though  the  momentum  algorithm  can  easily  be  used 
to  find  the  extreme  vertices  of  a  RW  momentum 
envelope,  identification  of  connected  points  must 
be  done  manually.  Automating  this  task  would 
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provide  an  extremely  rapid  creation  of  the  edges 
of  the  geometry. 

•  The  generation  of  the  quasi-maximum  space  for  a 

hybrid  array  currently  depends  on  a  very 
simplistic  method  examining  every  CMG  point, 
finding  the  nearest  RW  point  and  adding  the 
momentum  values,  subject  to  an  angle-error 

constraint.  Due  to  the  highly  irregular 

distribution  of  the  RW  points  and  the  very 

regular  distribution  of  CMG  points,  there  can 
often  be  a  miss  where  a  CMG  point  has  no 

corresponding  RW  point.  Due  to  the  geometrically 
simple  shape  of  RW  momentum  spaces,  it  should  be 
possible  to  determine  the  momentum  value  at  any 

point  on  the  surface  and  add  that  value  to  the 

CMG  point. 

•  Further  analysis  should  be  conducted  to  determine 
the  reasonableness  of  including  small  RWs  with  a 
CMG  array  to  allow  torque-free  elliptic 
singularity  avoidance  within  the  entire  CMG 
momentum  space . 

2 .  Torque  Analysis 

•  The  box/hypercube  singularity  free  search 
algorithm  misses  a  significant  amount  of  adjacent 
singularity  free  space  due  to  the  simplicity  of 
implementation  for  both  the  search  algorithm  and 
the  resulting  singularity  free  momentum  space 
analysis.  This  unnecessarily  restricts  CMG 
capability.  An  improved  search  algorithm  would 
restore  some  of  the  lost  performance. 

•  Due  to  the  numerical  nature  of  the  singularity 
free  space  created  by  the  state  space  analysis, 
another  method  must  be  developed  to  verify  that 
the  space  is  actually  singularity  free.  Ideally, 
this  could  be  done  analytically  as  it  would 
possibly  allow  an  analytic  determination  of 
allowable  gimbal  angle  ranges,  obviating  the  need 
for  a  numerical  analysis. 

•  Adding  limited  gimbal  angles  to  the  momentum 
inversion  algorithm  would  allow  the  singularity- 
free  momentum  space  to  be  created  with  the  same 
method  as  the  normal  maximal  surface. 
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•  Further  development  of  the  control  law  in 

Equation  (94)  would  allow  full  utilization  of  the 
singularity  free  space  provided  by  the  gimbal 
angle  limitations. 

3.  Reaction  Analysis 

•  The  value  of  the  variables  and  parameters 

involved  in  the  force  analysis  can  vary  across 
many  orders  of  magnitude.  To  prevent  the 
insertion  of  numerical  errors  when  performing 

arithmetic  between  very  large  and  very  small 

numbers,  scaling  should  be  introduced  to  the 
equations . 

•  The  simplified  model  may  not  accurately  model  the 

physical  MEDs .  Further  analysis  should  be 

performed  to  improve  the  simplifications. 

•  The  design  of  a  MED  array  could  be  optimized  for 

a  specific  mission  by  assuming  the  operating 
regime  of  a  spacecraft  and  designing  the  MED 

based  on  the  calculated  loads 
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